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Make speedier checks of recorders, control- 
lers and base or noble metal thermocouples 
in industrial plants with the new three-dial 
8686 Portable Millivolt Potentiometer. Fea- 
tures such as a central reading window... 
where measured values appear as a row of 
digits with a scale interpolation . .. simplify 
calibration of thermocouples and test meas- 
urements. The 8686 Potentiometer has: a 
wide operating range of —10.0 to +100.1 
mv and +1010 to +1020 mv for standard 
cell calibration; and a high accuracy of 
+(0.05% of reading +3uv) without refer- 
ence junction compensation, +(0.05% of 
reading + 6 nv) with ref. jet. comp. Write 


8686 Millivolt 
Potentiometer 


8690 Millivolt 


for Data Sheet E-33(1A). 


8692 and 8693 
Temperature 
Potentiometers 


If you want to make a variety 
of temperature measurments 
quickly with one flexible 
instrument, investigate the 
new time-saving 8692 Single- 
Range or 8698 Double-Range 
Temperature Potentiometers. 
Available in any of 24 inter- 
changeable temperature and 
millivolt ranges, these instru- 
ments read directly in degrees 
F or C on a scale 27%” long. 
Convenience features include: 
simplified range changes . . . 
only a screwdriver is needed 
to change a circuit panel, 
scale and binding post studs; 
automatic reference junction 
compensation . . . reference 
coil, built into circuit panel, 
compensates for thermocouple 
being used; accuracy . ‘ 
+0.2% of range. Write for 
Data Sheet ND42-33(1A). 
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Need a fast-operating, high-sensitivity, 
high-quality null indicator for use in re- 
search, testing and production checking? 
Here’s a new 9834 Guarded D-C Null Detec- 
tor having a short period of less than two 
seconds for source resistances up to 1000 
ohms, increasing to 4 seconds at 100,000 
ohms . . . ideal for measurements with 
guarded or unguarded potentiometers and 
bridges. Of rugged construction, this port- 
able, line-operated detector provides nu- 
merous convenience features which include 
four degrees of sensitivity, with a basic 
sensitivity of 0.2 uv/mm (0.3 uv/scale div.), 
and a noise level of less than +0.1 pv. 
Write for Data Sheet ED7(2). 
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Potentiometer 


Fast, accurate in-plant check- 
ing of thermocouples, record- 
ers and controllers is expe- 
dited by the new 8690 Port- 
able Millivolt Potentiometer. 
Design features such as “In- 
line Readout”, which permit 
measured values to be read in 
digits with a scale interpola- 
tion, accelerate temperature 
measurements. The 8690 pro- 
vides a wide range of —11.0 to 
+101.0 mv which eliminates 
input reversals; an accuracy 
of +(0.05% of reading + 20 
uv) without reference junc- 
tion compensation, + (0.05% 
of reading + 40 uv) with 
ref. ject. comp. Write for Data 
Sheet E-33(1B). 


9834 Guarded 
D-C Null Detector 


4992 Stenton Ave., Te Philadelphia 44, Pa. 
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A COMPLETE PRINTING SERVICE 


Goop PRINTING does not just happen; it is the 
result of careful planning. The knowledge 
of our craftsmen, who for many years have 
been handling details of composition, print- 
ing and binding, is at your disposal. For 
over seventy-five years we have been printers 
of scientific and technical journals, books, 
_ painters or . theses, dissertations and works in foreign 


languages. Consult us about your next job. 


LANCASTER PRESS, Inc. 
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Announcing publication of 


The Proceedings of the Symposium on 
—GAS-COOLED REACTORS- 


held February 10-11, 1960 


sponsored jointly by The Franklin Institute and the Delaware 
Valley Section of the American Nuclear Society. 


Twenty-one papers are featured by well known 
U. S. research people and several foreign scientists 
who have reviewed the developments in this new type 
of atomic power plant. The symposium aimed to 
clarify the problems which must be solved before gas 
cooling can be used to its fullest potential. 


Sir William Cook, C.B., of the U. K. Atomic Energy 
Authority reviewed the British program in gas-cooled 
reactors; a representative from the Swiss firm, Brown, 
Boveri & Cie., discussed their work on a pebble bed 
reactor which they are building for the West German 
Government. Representatives of the Australian Atomic 
Energy Commission and from West Germany pre- 
sented programs being carried out in their respective 
countries. Discussions are included at the end of 
most papers. 


The sessions were devoted to the following topics: 
I Status of Major Gas-Cooled Reactor Programs 


I! Fuel Elements and Materials 
111 Components and Loop Studies 


Issued as Monograph No. 7 under the auspices of the JOURNAL OF THE FRANK- 
td mahi the Proceedings are now available. Send your order (at $5.00 per 
copy) to 


JOURNAL OF THE FRANKLIN INSTITUTE 


20th & Parkway 
Philadelphia 3, Pa. 
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HAROLD S. BLACK, 
LAMME MEDALIST 


JEWS HOPEFULLY 
THE HITHERTO 
UNATTAINABLE 


A MAN WINS A MEDAL...AND 
STRENGTHENS A PHILOSOPHY 


The search for the “hitherto unat- 
tainable” sometimes ends in strange 
places. 

For years Bell Laboratories engi- 
neer Harold S. Black pondered a 
problem: how to rid amplifiers of the 
distortion which unhappily accumu- 
lated as signal-transmission paths 
were made longer and amplifiers were 
added. There had been many ap- 


proaches but all had failed to provide 
a practical answer. 


Then one day in 1927 the answer 
came—not in a research laboratory, 
but as he traveled to work on the 
Lackawanna Ferry. On a newspaper, 
Mr. Black jotted down those first ex- 
citing calculations. 


Years later, his negative feedback 
principle had revolutionized the art 
of signal amplification. It is a prin- 
cipal reason why telephone and TV 
networks can now blanket the coun- 
try, the transoceanic cable is a reality, 
and military radar and missile-con- 
trol systems are models of precision. 


For this pioneer achievement, and 


for numerous other contributions to 
communications since then (some 60 
U. S. patents are already credited to 
him), Mr. Black received the 1957 
Lamme Medal from the American In- 
stitute of Electrical Engineers. He 
proved that the seemingly “unattain- 
able” often can be achieved, and thus 
strengthened a philosophy shared by 
all true researchers. 


He is one of many Bell Telephone 
Laboratories scientists and engineers 
who have felt the challenge of teleph- 
ony and have risen to it, ranging 
deeply into science and technology. 
Numerous medals and awards have 
thus been won. Two have been Nobel 
Prizes, a distinction without equal in 
any other industrial concern. 


Much remains to be done. To cre- 
ate the communication systems of the 
future, we must probe deeper still for 
new knowledge of Nature’s laws. We 
must continue to develop new tech- 
niques in switching, transmission and 
instrumentation for every kind of in- 
formation-bearing signal. As never 
before, communications offer an in- 
spiring challenge to creative men. 


BELL TELEPHONE LABORATORIES 


WORLD CENTER OF COMMUNICATIONS RESEARCH AND DEVELOPMENT 
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OUR UNBALANCED BIAD * 


BY 


CHAUNCEY D. LEAKE! 


The title of my talk, “Our Unbalanced Biad,” contains a word not 
to be found in our dictionaries. You may, nevertheless, realize quickly 
what I am to talk about, when I suggest that there are two arms to our 
intellectual balance. Part of our continuing job of learning is to keep 
these long arms of intellectual effort in balance 
one hand and the sciences on the other. 

In our time we are condensing the long respected Greek triad of 
logics, ethics and aesthetics into the biad of the sciences and the humani- 
ties. Are we wise? Is the balancing of our scientific disciplines with 
our humanistic studies to be enough? 


THE GREEK TRIAD 


One may find existential relations between the concepts of logics, 
ethics and aesthetics. In order to gain satisfactions in living, it is 
necessary first to find somewhat of what we call the ‘‘truth”’ about our- 
selves and the universe surrounding us. This is the business of logics. 

In the light of such knowledge as we acquire, we are then prepared 
to examine the motives which determine our conduct and our relations 
with other people, in seeking the goals we prefer in our hierarchy of 
values. Here is the problem for ethics. 

Then comes the question of how best to apply our knowledge to the 
achievement of our purposes, to reach the satisfactions we crave. This 
is the function of aesthetics. 


* Address delivered at the Annual Meeting of The Franklin Institute on Wednesday, 
January 20, 1960. 

1 President, American Association for the Advancement of Science; Professor of Pharma- 
cology, The Ohio State University, Columbus, Ohio. 

(Note—The Franklin Institute is not responsible for the statements and opinions advanced by contributors in 
the JouRNAL.) 
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In our current,terminology, the logics have become our sciences, the 
aesthetics our arts. What has happened to the ethics? Have these 
evaporated into our ‘‘value judgments,” or are they condensed with 
our arts into the ‘humanities’? 

Many ways of obtaining truth have been tried, from oriental em- 
phasis on subjective experience, including mystical revelation, to occi- 
dental preoccupation with objective observation. Aristotle’s great 
achievement was to systematize principles of reasoning. His deductive 
formal logic became so powerful that it was the logic, so that the word 
itself was used only in the singular. Aristotelianism dominated occi- 
dental intellectual affairs for centuries. It was brilliantly exploited by 
the Christian Fathers, and became the foundation of scholasticism. 

Other forms of logic were not revived until recently. Inductive 
logic, frequently used in the natural sciences, grew with the Renaissance. 
A logic of science, with verifiability by measuring, gradually evolved. 
Symbolic logic has developed with mathematics. Other non-Aris- 
totelian logics are being explored. 

Interestingly, modern neurophysiology is confirming psychoana- 
lytical insight by finding analogical processes associated with coding 
mechanisms in the functioning of our brains. These seem to be prime 
in creativity. Analogy involves recognition of potential relationships 
between what has hitherto been conventionally thought to be unrelated. 
Free association, or free conscious analogy, may be the greatest stimulus 
to creativity in both the arts and sciences. Its honest use may itself 
free us from our stereotyped conventionalities. 

The great merit of the scientific method is that verification by refer- 
ence to agreed-upon standards of measurement establishes a basis for 
agreement on what constitutes the ‘“‘truth.’’ It is interesting how the 
psychological factor of voluntary agreement becomes paramount in 
science. Even the standards used for measuring have validity merely 
by agreement. 

It is this verifiability that assures the strength of modern science. 
However, in our naive awe at this enormously powerful tool, let us 
remember that wisdom is more than verifiable and agreed-upon knowl- 
edge of ourselves and our environment. There still remain the prob- 
lems of what to use this knowledge for, and how best to use it for the 
purposes we may agree upon. 

In approaching the ethical problem, the ancient Greeks saw the 
conflict between conduct directed toward self-gratification on the one 
hand, and the welfare of society on the other. The relation of indi- 
viduals to the societies in which they live remains vastly complicated 
with uncountable and ever-changing factors. Yet both old and new 
democracies depend on voluntary agreement among the people in a 
community as to appropriate goals toward which they might willingly 


work together. 
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The aesthetic problem involves judgment in choice of the most 
fitting or appropriate way by which our agreed-upon (or scientific) 
knowledge of ourselves and our environment can be applied to the 
accomplishment of our agreed-upon purposes. Judgment is sharpened 
only by practice or experience, as every artist or engineer or physician 
learns. <A successful aesthetic solution calls for discerning taste in the 
fitness or propriety of things or ideas. Voluntary social agreement is 
more difficult to obtain in aesthetics than in logics or ethics, and 
when it is reached, there tends to be conventionalization, with status 
stratification. 

The triadic relationship among logics, ethics and aesthetics is 
mediated by the neglected factor of voluntary agreement, and may be 
expressed in terms of objectives: the function of logics is to establish 
truth; the goal of ethics is to obtain good ; and the purpose of aesthetics 
is to achieve beauty. The true, the good, and the beautiful thus com- 
prise a trinity which retains its appeal and its power. There seems to 
be an indication that the truth itself points toward that which is good 
and beautiful. 

The Greek triad, relating logics, ethics and aesthetics, probably 
reached its clearest expression in the 5th Century B.C. under the in- 
fluence of Plato. Thinkers have played variations on the theme ever 
since. By no means have we fully explored its capabilities or its sig- 
nificance. Peculiarly, in our critical time, we are attempting to alter 
the stability of this triad, and perhaps foolishly to simplify it. 


THE CURRENT BIAD 


Slowly, imperceptively, almost unconsciously, our academic trend, 
which sets the direction of our intellectual life, has narrowed the 
traditional triad of logics, ethics and aesthetics to the current biad of 
the sciences and the humanities. Significantly, and perhaps to our woe, 
ethics seem either to have dropped out of social consideration, or to 
be somehow or other included among the humanities. 

Our age of science is a time of crisis. Our value judgments are in 
chaos. Old values are weak and new ones uncertain. Many of our 
cherished beliefs are being shaken by the rough demonstrable facts 
about ourselves and our universe. Yet it is part of our scientific spirit, 
as my revered teacher, Edwin Grant Conklin, emphasized, to be con- 
fident that unwanted truth is better than cherished error. Here is a 
supreme value judgment. 

Our occidental culture has withstood many crises. The ancient 
myths, now being so entertainingly examined by Robert Graves, may 
record the overthrow of matriarchal societies with the growth of patri- 
lineal tribes, as human males realized that they had something to do 
with continuing the race. Later legends indicate the decline of nature 
worship in the face of a vigorous Greek rationality, which in turn de- 
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teriorated under Roman luxury and administrative rigidity. Popula- 
tion pressures destroyed the materialistic stability of the Caesarean 
world, and into the darkness came a yearning faith in spiritual reward, 
which crystallized in an authoritarianism built on absolute ideals. 

The great Thomistic achievement in the 13th Century of a solid 
faith corresponding precisely with then demonstrable fact began to 
crack as the alleged facts were questioned by a skeptical humanism. 
Copernican theory, extended by Kepler, toppled the cosmology of 
antiquity; Harvey overthrew the neat Galenical system of how our 
animal bodies work; Galileo, Boyle and Lavoisier shook intellectual 
complacency with the verifiable measurements that started modern 
physics and chemistry; Ludwig and Helmholtz issued a manifesto of 
faith that all living processes are explainable in terms of physics and 
chemistry, with lasting reverberations in the Russia that became the 
testing ground for another manifesto; and Darwin roughly disturbed 
the comfortable religious beliefs of our great-grandparents a century ago. 

It is interesting to note how Newton’s rigid system of measured 
checks and balances, which even set a pattern for political and economic 
doctrine, as reflected in the United States Constitution, should yield 
first to the concept of gradual change through evolution as introduced 
by Darwin, and then to the more general idea of relativity as developed 
by Einstein. It is ironic that now we seem to have come to the cer- 
tainty of uncertainty ! 

As we pursue our probing of the world around us, we go from space- 
time confusion in the tiniest bits of organized energy in nuclear par- 
ticles to time-space bewilderment in receding nebulae and an unbeliev- 
able expanding universe, with infinity compounded by indefiniteness. 
We find no answer to the riddle of being. 

Life itself baffles us. We observe in the macromolecules of genes 
and viruses no clear dividing line between what lives and what doesn’t. 
We learn, following Virchow’s idea of studying disease from cells to 
societies, that organizational levels of living material involve factors 
not capable of physical or chemical explanation. Further we discover 
that we ourselves, in the very act of observing, modify what we ob- 
serve. The psychology of gestalt further breaks down the materialistic 
mansion we have built and come to fear. 

This then is another time of intellectual crisis. It is brought on by 
the very success of our scientific logic. Planck’s quantum theory, with 
the thermodynamics of Gibbs and Maxwell, and Einstein’s mathe- 
matical logic, have given us the power of nuclear energy, with a glimpse 
into the energy factors which may control life itself. The failure of our 
ethics appears in our distressing dilemma of decision: shall we use 
nuclear power for arrogant threats in daring mutual destruction, or 
shall we use it for the peaceful benefit of all people? 

Our social schizophrenia can be diagnosed from our aggressive prac- 
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tice in the face of our professed peaceful intent. Effective therapy 
directed toward bringing into balance the applications of our sciences 
with the implications of our humanities may depend on sharpening the 
fulcrum of our self-knowledge. Perhaps Socratic echoes may still 
stir us. 

Our dissatisfactions grow from our anxieties that our generous good- 
willed ideals of behavior cannot be realized in the crass materialism of a 
machine-like world. Such anxieties yield to frustrations as we begin 
to see that we may not get the sort of satisfactions we want. Then our 
frustrations focus in resentments against those we think are blocking 
our paths. Thus our scientists are blamed for our manifest distress, 
and an anti-intellectualism arises. Artistic endeavor becomes un- 
intelligible to the uninitiated ; blasé beatniks with their ladies in Levis 
squawk rebellion ; and conventional wisdom calls for more science in the 
name of security. 

Our international and individual crises are produced both by insufh- 
cient knowledge, especially of ourselves, and by inability through fear 
to join voluntarily in social agreements as to common purposes. A 
psychology of determent emerges, perhaps to stabilize a nobler psy- 
chology of responsibility. 


OUR LACK OF BALANCE 


We are thus beginning to sense our difficulties and our failures. In 
seeking by scientific means to obtain the truth about ourselves and our 
environment, we seem to have put too much emphasis on our environ- 
ment to the neglect of a merciless study of ourselves. In the relation- 
ship, then, between ourselves and our environment, we are out of 
balance. 

Lacking directional ethics, we cannot apply our scientific knowledge 
to purposes we can agree upon as worthy of universal social endeavor. 
Our science seems to become applied chiefly to making individual and 
national status symbols. What gadgetry is directed toward cleanliness, 
comfort and convenience, is well diluted by styling for false-fronting, 
and the waste is monumental. What profit longer life from medical 
science if there is so little chance yet for living it with graciousness, 
dignity, self-confidence and self-reliance? What social value has our 
science if it is applied so successfully to the possibility of destroying 
us all? 

The social impact of science has been tremendous. We have brought 
many frightening problems for which we are quite unprepared. Of 
universal concern are such as: (1) radiation hazards; (2) over-popula- 
tion dangers; (3) harmful pollution of airs and waters; (4) barriers of 
misunderstanding and distrust resulting from inequalities in the dis- 
semination and appreciation of scientific knowledge; (5) fears of a 
stereotyped scientific authoritarianism; and (6) industrial uglification. 
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Many of the applications of our scientific knowledge thus bring 
further imbalance in our cultural endeavor. If we are to build a 
balanced biad for measuring the success of our living adventure, we 
must emphasize the humanities along with the sciences. Our rapidly 
growing verifiable knowledge about ourselves and our environment 
must be balanced by an equal growth in wise judgments in the applica- 
tion of this knowledge to the attainment of the good and the beautiful. 

Yet even this balance is not enough, for we seem to be drifting with- 
out agreed-upon purposes. ‘‘All we like sheep have gone astray,” and 
instead of recognizing our individual responsibilities, we wait for some- 
one else to tell us what to do. Like waiting for Godot, we wait for 
God, not realizing that within us always is our God, who unites us to 
all that is not us, to all the rest of the universe, of which we are all a part. 

Maybe the attempt to build a balanced biad of the sciences and the 
humanities is not enough. Ethics is lacking to stabilize the balance. 
At best a balanced biad is in unstable equilibrium: it can so easily get 
out of balance by a shift in interest either on the scientific or the human- 
istic side. We crave stability. Perhaps some of the traditions of our 
culture, such as those nourished at The Franklin Institute, can help in 
lighting our way. Perhaps we might well give up the idea of maintain- 
ing a balanced biad between the sciences and the humanities, and return 
to the more stable and fundamental Greek concept of the triad of logics, 
ethics and aesthetics. 

This would be to our advantage. We now have a solid basis for 
logics in the verifiability of science. We have a keen appreciation of 
aesthetics in humanistic endeavor. We need, however, to reach agree- 
ments on appropriate social goals. These may be supplied by ethical 
studies. 

It is regrettable that philosophy and ethics have generally de- 
teriorated in our academic environments during the last half century. 
This has coincided with the rise of the materialistic applications of our 
science. We might greatly benefit were philosophical considerations 
and discussions to come more widely into our daily lives and activities. 
It is here especially that our liberal arts colleges can show the way, even 
for our huge state universities. ' 

Perhaps the philosophers are partly at fault. They seem to have 
developed, quite as scientists and artists have done, a peculiar jargon 
and vocabulary of their own, which is not understood or comprehended 
by people generally. Science, art and philosophy cannot flourish in 
secret. If they are to have meaning for humanity, they must be part of 
humanity. 

THE HOPE FROM UNDERSTANDING 

Often in the challenge of an intellectual crisis, despair may take the 
place of hope and faith. Understanding, however, as Leonardo da 
Vinci proclaimed, is the basis for love, and in such understanding there 
may come renewed hope and faith. 
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As our world-wide communications continue to improve, we may be 
sure that the traditional aspirations of people everywhere for a better 
life will push for the balanced triad of logics, ethics and aesthetics: our 
science will become more widely understood, and as our inherent love of 
beauty exerts itself, our applications of scientific knowledge may be ever 
more wisely directed toward social welfare and the good of people. 
We are witnessing a new resurgence of enthusiasm for understanding 
science, for achieving worthy purposes in human welfare, and in pre- 
serving the natural beauties of our world, even in our increasing 
industrialization. 
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COMPARISON OF YIELD CONDITIONS FOR CIRCULAR 
CYLINDRICAL SHELLS * 


BY 
ANTONI SAWCZUK'! AND PHILIP G. HODGE, JR.’ 


ABSTRACT 
The yield condition for a rigid-plastic circular cylindrical shell subjected to 
axially symmetric loading is described in terms of the circumferential hoop stress n 
and the axial bending moment m. Different yield conditions are obtained depending 
on whether the shell is of homogeneous or ideal sandwich construction, and on whether 
the shell material satisfies Tresca’s or Mises’ yield condition. Previous work has been 


limited to Tresca’s yield condition. 

The example of a long shell subjected to a ring of pressure at its midpoint is con- 
sidered for each of the four shell yield conditions as well as for a simple “limited inter- 
action” approximation. Complete stress and velocity fields are determined, the 
limit load is found, and the extent of the plastic region is determined. ‘The results 
are assessed in terms of the theorems of limit analysis. 


1. INTRODUCTION 


Previous analyses of plastic shell problems have usually started from 
the premise that the shell material satisfies Tresca’s yield condition of 
maximum shear stress. However, experimental evidence generally in- 
dicates that the actual yield condition is a very complex one and that 
the Mises yield condition of maximum octahedral shearing stress is a 
somewhat better approximation to it. Therefore, it seems of some 
interest to investigate a simple shell problem under the hypotheses of 
various yield conditions and to compare the results. 

In addition to considering different material yield conditions, it is 
also of interest to consider different types of shell construction. Prob- 
ably the most common type of shell is one of constant thickness with 
uniform properties across the thickness. However, it is sometimes 
mathematically convenient to consider an idealized sandwich shell 
consisting of two ideally thin sheets separated by a core of finite 
thickness. 

The particular problem chosen for consideration is that of a very 
long circular cylindrical shell, subjected to a single ring load. The shell 
is assumed to be sufficiently long so that the support conditions at the 
end are immaterial. In this case, the plastic action is confined to a 
finite region symmetric about the ring of force. The complete problem 


* The results presented in this paper were obtained in the course of an investigation 
sponsored by the Office of Naval Research. 

! Post-Doctoral Fellow, Illinois Institute of Technology, on leave from the Polish Academy 
of Sciences, Chicago, III. 

2 Professor of Mechanics, Illinois Institute of Technology, Chicago, III. 
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is to determine the extent of this plastic region, the plastic collapse load 
(the magnitude of the force necessary to cause plastic action), and the 
stress and velocity distributions throughout the plastic region. 

A complete solution of this problem has previously been given (1)? 
for a simplified limited interaction curve‘ which has no direct physical 
meaning. In addition, the plastic collapse load has been computed for 
the Tresca yield condition, for both the uniform and sandwich shells (1). 
The present paper derives the complete stress solutions for the four 
combinations of yield conditions and shell structure and the velocity 
solutions for the sandwich shells. 


Rigid region Plastic region Rigid region 


2H 


(d) 


Fic. 1. Cylindrical shell with ring load. 


(a) Loaded shell 

(b) Stress resultants 

(c) Section of uniform shell 
(d) Section of sandwich shell 


Section 2 of this paper gives the basic equations of the problem and 
the interaction curves for the various theories. A general method of 
solution is derived in Sec. 3. Section 4 then presents formulas for the 
various shells considered. Finally, the different solutions are com- 
pared in Sec. 5. 


3 The boldface numbers in parentheses refer to the references appended to this paper. 
4 So-called because the restrictions on direct force and moment are assumed to be inde- 
pendent of each other. 
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2. BASIC EQUATIONS 


The basic equations of cylindrical shell theory are well known (see, 
for example, (2)). The problem is defined in terms of the radial velocity 
W and three stress resultants per unit length: the circumferential direct 
stress No, the axial moment M,, and the shear force Q. Sign conven- 
tions for the stresses are defined in Fig. 1; W is positive outward. We 
find it convenient to use the following dimensionless quantities : 


= Ni/No m= M,./M, w= W/A 
x= Mo) 2X 
My = uniform shell 
sandwich shell 


(1) 


| 


-06 -04 -02 16] 02 04 


Fic. 2. Interaction curves. 
——————— Uniform Mises shell 
——- ——-_ Sandwich Mises shell 
—-— Uniform Tresca shell 

Sandwich Tresca shell 

——- - —— Limited interaction curve 


The coordinate is measured from the point of application of the ring 
load; primes will be used to represent derivatives with respect to x; 
a» and o,* are the yield stresses of the uniform shell and sandwich shell 
materials, respectively ; shell dimensions are shown in Fig. 1. 

The behavior of a rigid perfectly plastic cylindrical shell is defined 
by two equilibrium equations, an interaction curve, and a flow law. 
In dimensionless terms and with no distributed surface loads present, 
the equilibrium equations are 


qgt+n=0 (2) 
= (3) 
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It is assumed that the effect of the shear force on plastic yielding is 
negligible. Therefore, a description of the plastic behavior can be 
given in terms of m and m only. It is convenient to give this descrip- 
tion in a stress-resultant-space with coordinates m and m. 

We begin with the interaction curve which plays the same role in 
stress-resultant-space as the yield curve in stress-space. Indeed, 
beginning with a given shell construction (uniform or ideal sandwich) 
and a given material yield condition (Mises or Tresca), one can derive 
the corresponding interaction curve in terms of mand m. Alternatively, 
one may assume the form of interaction curve directly as a property of 
the integrated shell based, for example, on experimental results or 
mathematical convenience. Some specific curves are shown in Table I 


TABLE |[.—Jnteraction Curves. 


Shell 

Structure Condition Ref. Interaction Curve 

Uniform Mises (3) n = + tan p log tan (p/2) 

m = + (2/v3)[tan® p log tan (p/2) + sec p] 
O<p<cr 


sin p m = — (2/v3) cos p 
O<p<2r 


Sandwich Mises (3, 4) n 


ll 


Uniform Tresca (1, 2) + 4n(1 + n) - 


Sandwich Tresca (5, 2) 


Limited (t,.2) n 
Interaction 


and Fig. 2. In every case the interaction curve is closed, convex, and 
symmetric with respect to both m and n. 

If the stress point for a particular value of x is within the interaction 
curve, no plastic flow is possible ; thus for our ideal rigid-plastic material 
the shell is rigid. 

If the stress point is on the interaction curve, plastic flow may occur 
according to the plastic potential flow law. This flow law depends not 
only upon the state of stress but also upon the assumed interaction 
curve. Let e, and k, be the dimensionless unit circumferential exten- 
sion rate and axial curvature rate respectively, of the middle surface of 
the shell and superpose axes k,, és on the m, m axes in Fig. 2._ The flow 
law then states that the ‘‘strain-rate vector’’ with components k,, é 


’ 
: 

1 <n < — 1/2 
m=+1 —1/2<n<1/2 
oe m = + 4n(1 — n) 1/2<s2<1 
m=2+1 —1/2<n<1/2 
m= + 2(1 — n) 1/2<n<1 4 
=+1 ~ises1 
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must be directed along the outward normal to the interaction curve at 
the stress point in question. In terms of the present variables, 


Therefore, if dm and dn represent increments along the interaction curve 
the plastic flow rule states that 


w''dm + wdn = 0. (4) 


Finally, since the material does not exhibit strain hardening, the 
stress point can never lie outside of the interaction curve. 

In general, the stress distribution in a rigid-plastic material may not 
be uniquely determined, but the plastic collapse load is necessarily 
unique. It follows that the statement of the basic boundary value 
problem must be carefully made. We are asked to determine four 
functions m, n, g, and w and a load magnitude go such that the following 
requirements are all satisfied : 


(a) Eqs. 2 and 3 are everywhere valid ; 

(b) the stress point (m, 2) is everywhere on (plastic) or inside of 
(rigid) the interaction curve ; 

(c) in a plastic region, Eq. 4 is valid; 

(d) in a plastic region, the strain rate vector (w’’, w) is directed out- 
wards from the interaction curve; 

(e) in a rigid region, w = 0; 

(f) all stress and velocity boundary conditions are satisfied ; 

(g) the moment m, shear force g, and velocity w are everywhere 
continuous; 

(h) the slope w’ is continuous except that if m takes on its extreme 
possible value, w’ may have a discontinuity (hinge circle) of the 


same sign as m. 
A physical consideration of the problem suggests that the plastic 


deformation will be confined to a central portion —x, < x < x, of the 
shell and that > Oin this region. Therefore, the boundary conditions 


for the plastic region are 
q = qo (5) 
w' = 0 minimum at x (6) 
= maximum at x (7) 


(8) 
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3. GENERAL SOLUTION 


To obtain a general solution to the problem defined in Sec. 2, we 
express all quantities in terms of a parameter p chosen so that dm/dp >0. 
In particular, m and n are known functions of p given by the parametric 
representation of the interaction curve: 


m=m(p) n=n(p). 

Equations 2 and 3 then become 
(dq/dp)p’ +n = 0 (10a) 
(dm/dp)p’ = (100) 


Elimination of p’ from Eqs. 10a and 104 then leads to an equation whose 
solution can be written 


n(dm/dp)dp. (11) 


(9) 


Here pp is the value of p at x = 0, and the boundary condition (Eq. 5) 
has been incorporated in the solution. 

If g, and p, denote the values of g and p at x = x;, Eq. 11 can be 
solved for the collapse load qo: 


f° n(dm/dp)dp. (12) 


The lower bound theorem of limit analysis (2) States that any stress 
field which satisfies conditions (a) and (d) and the stress portions of (f) 
and (g) will provide a lower bound on the collapse load. Assuming, for 
the moment, that the stress field can be continued into the rigid region, 
any choice of po, 1 and q; will provide a lower bound on qo. Holding 
q: fixed and noting that the integrand in Eq. 12 is positive, we see that 
the highest lower bound is obtained when p; is a maximum and fo a 
minimum. This implies that m isa minimum at x = 0 and a maximum 
at x = x, But if m is a maximum at x, and m’ is continuous there, 
=q=4q: = 0. Therefore, Eqs. 12 and 11 become 


E f (dm/dpyap |" (13a) 


PL 1/2 
E f (dm /ap)dp| (13) 
With g known as a function of p from Eq. 13), Eq. 10 can be written 


(dm/dp)(dp/dx) = q 
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and solved explicity for x: 


x= f’ (1/q)(dm/dp)dp. 


In particular, the plastic-rigid boundary is given by 


f (1/q)(dm/dp)dp. (13d) 


In general, Eq. 13¢ must be integrated numerically. Once this has 
been done, m, n, and g may all be computed as functions of the coordi- 
nate x from Eqs. 9 and 130. 

If the shell is sufficiently long, it is evident that some solution of 
Eqs. 2 and 3 for x > x; can be found such that the prescribed support 
conditions at the end of the shell will be satisfied and the stress point 
will always be on or inside of the interaction curve.’ Therefore, the 
proposed solution satisfies all requirements on stresses. 

Next, we shall show that a velocity field of the form 


(14) 


satisfies all requirements on velocity. Indeed, it follows from Eqs. 2 
and 3 that 


q''(dm/dp) = — n'(dm/dp) = — (dn/dp)p'(dm/dp) 
= — m'(dn/dp) = — q(dn/dp) 


so that Eq. 14 satisfies Eq. 4. Further, m is such that the boundary 
conditions (Eqs. 6 and 7) are satisfied, and Eq. 14 satisfies the boundary 
condition (Eq. 8). Finally, it follows from Eq. 13d that g, and hence 
w, is always positive; therefore the strain-rate vector is directed out- 
wards from the top side ” > 0 of the interaction curve. 

Therefore, the proposed solution satisfies all of conditions (a-h) 
and hence our hypothesis is verified. 


(15) 


4. EXAMPLES 


We consider first the sandwich shell with Mises’ yield condition. 
The interaction curve in Table I is already in the proper parametric 
form, and 


 pi=t. (16) 

Substituting from Table I and Eqs. 16 into Eqs. 13 we obtain 
go = (2x/ = 1.905 (17a) 
q = (1/v3)[2(4 — p) + sin 2p}! (176) 


5 The question of how long the shell must be and the modification necessary for short shells 
has been dealt with for a particular interaction curve (2, 6) and will not be considered here. 


( 
33 
| 


May, 1960.] YIELD CONDITIONS FOR CIRCULAR SHELLS 369 


= f sin p[2(x — p) + sin (17c) 


x, = 3.467. (17d) 


For the uniform shell with Mises’ yield condition, the appropriate 
portion of the interaction curve is given by the lower sign in Table I. 
Equations 16 are still applicable, so that Eqs. 13 become 


Go = 21/2.3-3/4.¢ = 1,949 (18a) 


g = — (2/n*){2 tan* p[log tan (p/2) 
+ tan psec p log tan (p/2) — tan p 
sec p log tan (p/2)dp} (18d) 
__ _ 4 p[log tan (p/2) + cos p] 
x= q(p)cos? p dp. (1&¢) 


The numerical value of x; was not calculated for this case. 

The relevant portion of the interaction curve for the uniform shell 
with Tresca’s yield condition consists of the two arcs BEF, FGK in 
Fig. 2. These may be written in parametric form: 


n=1+p m-=4p(1 +?) 
n=1-p m= 4p(1—?) 


Therefore, Eqs. 13 lead to 


—1/2<p<0 
0<p<i/2 


(19) 


go= (10/3)'/# = 1.826 (20a) 


q = (1/v3) (5 —24p — 36p? — 16p*)'” —1/2<p<0 (20) 
q= (1/v3)(1—2p) (5 —4p)!" 0<p<i/2 

—1/2 (20c) 
x =4v3[0.1818+ V5 — (5—4p)!/7] O<p<t/2 
x1 =2.379. (20d) 


For the sandwich shell with Tresca’s yield condition, the pertinent 
part of the interaction curve is the broken line BFK in Fig. 2. This 
maybe written parametrically in the form 


n=1+p -—-1/2<p<0 
n=1-p m=2p 0<p <1/2 


(21) 
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Substitution of Eqs. 21 into Eqs. 13 then furnishes the results 


(1/v2)(3 — 8p — 
(1 v2) (3 Sp 


v2 sin=' [2(p + 


x = 2.532 — v2 cosh! 2(1 — p) 


x; = 3 


5 


Finally, for the limited interaction curve, the relevant portion 
can be written 


= m= p —-1<p<l (23) 


so that Eqs. 13 become 
(24a) 
(24b) 
(24c) 
x, = 2. (24d) 


The complete results for the limited interaction curve and some of 
the results for the Tresca yield condition were previously obtained by 
Drucker (1). 


—Comparison of Collapse Loads. 


1 2 3 4 5 


TABLE II. 


Difference Difference, % 


from 
Shell Yield Uniform Lower Upper 
Structure Condition qo Mises Actual Bound Bound 
Uniform Mises 1.949 0 0 0 0 
Sandwich Mises 1.905 — 0.044 —2.3 —2.3 2.6 
Uniform Tresca 1.826 —0.123 —6.3 —6.3 7.9 
Sandwich Tresca L732 —0.217 —11.1 —11.1 
Limited Interaction 2.000 0.051 2.6 —20.0 18.2 
Elastic 
(v = 0.3) 1.289 — 0.660 — 33.8 — — 


5. COMPARISON OF RESULTS 


We discuss first the collapse loads predicted by the various theories. 
For the sake of argument, we assume that the uniform shell with Mises’ 


‘ 
(J. F. 1. 
eat: go = V3 = 1.732 (22a) ae 
(226) 
q = 0<p<1/2 
(22c) 
a O<p<1/2 
(22d) 
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yield condition represents the actual problem to be solved, and that all 
other theories are to be regarded as approximate ones. The first 
column of Table II collects the collapse loads from Sec. 4, the second 
column shows the differences from the actual collapse load, and column 
three shows the percentage difference. 

In view of the lower bound theorem of limit analysis (2), if the 
collapse load is computed according to two interaction curves, one of 
which lies wholly interior to the other, the collapse load according to 
the interior curve is necessarily smaller than the other. It follows 
immediately from Fig. 2 that the collapse load for any problem com- 
puted for the sandwich Mises, uniform Tresca, or sandwich Tresca shell 
is a lower bound on the actual collapse load. 


Sandwich Uniform 
Mises +04 Tresca 


06 04 02 


02 04 06 O8 


Limited Sandwich 
Interaction o4+ Tresca 


Fic. 3. Approximate interaction curves. 


If the interaction curve is uniformly expanded by any factor, the 
corresponding collapse load will be multiplied by the same factor since 
the dimensioned collapse load is directly proportional to the yield stress. 
If, for example, we expand the interaction curve for the uniform Tresca 
shell by a factor of 2/v3 = 1.155, the new collapse load will be (1.155) 
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(1.826) = 2.108. However, as shown in the first quadrant of Fig. 3, 
such an expansion results in a curve wholly exterior to the actual one. 
Therefore, 2.108 is an upper bound on the actua! collapse load. In 
other words, by considering only the uniform Tresca shell, we can obtain 
both lower and upper bounds on the uniform Mises shell. 

vai Similar reasoning can be applied to each of the other approximate 
interaction curves. The various quadrants of Fig. 3 show the other 
curves scaled to be interior and exterior to the curve for the uniform 
Mises shell. The resulting bounds, shown as percentage differences 
from the actual collapse load of the uniform Mises shell, are given in 


Fic. 4. Moment distribution. 
— Sandwich Mises 
Uniform Tresca 
Sandwich Tresca 
—— - —— Limited interaction 
————— Elastic 


columns 4 and 5 of Table II. It is interesting to note that whereas the 
bounds on the error for the limited interaction curve are quite large, for 
this particular example the error is very small. Considering the mathe- 
matical simplicity of the solution for the limited interaction curve, it 
may frequently be worth using to obtain approximate solutions. 

Comparisons of the stress and velocity solutions are made in Figs. 4, 
Sand 6. In each case, the quantities are plotted at the corresponding 
collapse load, 
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Fic. 5. Hoop stress distribution. 


—— —— Sandwich Mises 
— — — Uniform Tresca 
Sandwich Tresca 
Limited interaction 
Elastic 


Figure 4 shows that there is relatively little difference between the 
bending moment distributions in the various cases. Although the 


Tresca condition theories predict a smaller plastic region, the moments 
within that plastic region are very close to the Mises sandwich shell. 


Fic. 6. Shear and velocity distribution. 


—— —— Sandwich Mises 
— Uniform Tresca ; 
* Sandwich Tresca 
— Limited interaction 
—————— Elastic (shear) 
———-— Elastic (velocity) 
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It appears reasonable to expect the actual solution to be similar to the 
approximate ones. 

Figure 5 for the hoop stress tells quite a different story. This is to be 
expected in view of the fact that at the end points of the plastic regime 
n = 0, 1/2, or 1 according to the Mises, Tresca, and limited interaction 
theories, respectively. Indeed, considering the differences in the bound- 
ary values, the values of ” for the various theories are surprisingly close 
together throughout most of the plastic range. Here again, we would 
expect the actual distribution to lie very close to that for the sandwich 
Mises shell. 

Figure 6 shows the shear force and velocity distributions which are 
proportional according to Eq. 14. For the Mises shell the slope is zero 
at x = 0 and x = x,, whereas there are hinge circles of discontinuous 
slope according to the other theories. Despite this difference in end 
conditions, the resulting deformed shapes are very similar through most 
of the plastic region. 

For comparative purposes, the elastic solution at the maximum 
elastic load is also shown in Figs. 4, 5, and 6. It is qualitatively quite 
different from the plastic solutions. This difference indicates that a 
considerable stress redistribution takes place during the stage of con- 
tained plastic deformation. In Fig. 6, of course, the elastic curves for 
shear force and velocity are quite different from each other since Eq. 14 
does not hold in this case. 
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SATELLITE PERTURBATIONS FROM EXTRA-TERRESTRIAL 
GRAVITATION AND RADIATION PRESSURE 


BY 
F. T. GEYLING' 


ABSTRACT 


The equations of motion of a satellite are written in terms of displacement com- 
ponents relative to the unperturbed, elliptic orbit. A moving system of coordinates 
is used which consists of an orthogonal triad whose origin is always located at the 
nominal satellite position on the elliptic orbit. This approach differs from the con- 
ventional one by way of Lagrange’s planetary equations in that the primary results 
are expressions for the satellite position rather than the osculating orbit elements, 
which seems an advantage for practical applications. Calculations are carried out 
for the effects of perturbations in the initial conditions of satellite motion and per- 
turbations due to oblateness, lunar and solar gravitation, and solar radiation pressure. 
In the numerical examples special attention has been devoted to the performance of a 
thin mylar balloon acting as a passive communication satellite. 


INTRODUCTION 


In the field of satellite perturbation analysis by far the largest 
amount of work has been based on Lagrange’s planetary equations 
where the equations of motion are reduced to six first order differential 
equations in the orbit parameters for which a first order approximation 
of the solution is obtained by the method of variation of parameters. 
The resulting expressions completely describe the time dependence of 
the osculating orbit and, in fact, of the actual orbit if one cares to carry 
out the necessary integrations which yield the perturbations in satellite 
position. The latter kind of information may be of chief interest to 
space technologists who classify themselves as amateurs by astronomers’ 
standards but who nevertheless need to know where a radar is to be 
pointed in order to track a vehicle whose motion has been characterized 
by previous measurements. In particular it should be convenient to 
know the extent to which the various physical perturbations cause a 
vehicle to depart from a position that can be forecast by simple elliptic 
orbit calculations for the unperturbed case. This idea of using an 
elliptic orbit as a reference is an obvious‘one and has been carried into a 
variety of integration schemes such as Hansen’s method, but to the 
present writer’s knowledge the existing literature on the motion of 
satellite vehicles contains very few analytic realizations of this approach. 

Specific calculations have been carried out for cases where the 
nominal, that is, unperturbed, orbit is circular and a generalization of 
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the results to elliptic nominal orbits will be given elsewhere.2 Since the 
effects of lunar and solar gravitation as well as radiation pressure are of 
some current interest in connection with communication satellites, many 
of the numerical examples have been tailored to this application. 


1, THE EQUATIONS OF MOTION 


For a dynamical system whose Lagrangian has been found, the 
Hamiltonian is defined as 


H(p,q,t) = Laipi-—L (1) 


where the generalized momenta are 


aL 
Og: 


By means of this last relation, the Hamiltonian can be written as a 
function of g, p, t. The canonical equations of Hamilton are then 


Ogi 


One advantage of the Hamiltonian approach to dynamic problems 
lies in the fact that the Hamiltonian may contain time explicitly, which 
can be due to a time-dependent coordinate system, time-dependent 
constraints, or a time-dependent potential function, V, in the La- 
grangian. A dynamical system whose /7 contains ¢ explicitly is said to 
be rheonomic. 

For the present work we find it attractive to select a moving coordi- 
nate system whose origin always coincides with the satellite position in 
unperturbed elliptic motion. This nominal satellite position is O’ in 
Fig. 1 where the plane of the ellipse is inclined by the angle i to the 
plane x, y. Any perturbations in the initial conditions or on the path 
will then result in satellite displacements about this point. A solution 
in this form is obviously convenient in determining satellite positions 
for tracking purposes. On the other hand one can still calculate secular 
perturbations for the osculating ellipse such as apsidal, anomalistic, 
and nodal precessions from these solutions without difficulty. We 


2 F. T. Geyling and H. R. Westerman, ‘“The Dynamics of Space Vehicles,”” Reading (Mass.), 
Addison-Wesley, in preparation. 
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describe the displacement from O’ in terms of u, v, w (Fig. 1) where the 
u-axis points radially outward from the ellipse, the v-axis at right 
angles thereto in the plane of the nominal orbit and in the direction of 
anomalistic motion, and the w-axis normal to that plane according to 
the right-hand rule. 
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In the x, y, z system of coordinates of Fig. 1 we have for the kinetic 
and potential energy of the moving body 


and 


GEm 
r 


V 


+7 
where 


m = mass of satellite 


E = mass of Earth 

G = universal gravitational constant 

r= [x2 + + 

V = whatever perturbative potential is being considered in addition 


to the first term in V, which of course represents a spherical 
Earth.’ 


Dots indicate time derivatives. 


§ A theorem of Poincare’s shows that for small perturbations the law of superposition holds, 
that is, various physical effects may be considered individually and their effects on satellite 
motion superimposed to describe the resultant orbit. 


: 
\ 


378 F. T. GeyLinc (J. F. I. 


According to the procedure which was just outlined, we first express 
7 and V in terms of u, v, w by the relations 


(u +r.) cos — 
—wsint + (vu sin @ cost + v cos 6 cost 
= wcost + (wu + sin @sinz + vcos sini 


where @ is the polar angle from the x-axis (which is the line of nodes in 
this case) to the nominal satellite position and r, the radial distance to 
that point on the nominal ellipse. We then obtain from (2) for the 
momenta 


= +r, — v6), p. = (6 + 7.6 + ub), Po = ti 


(where the right-hand sides are obviously the inertial velocity com- 
ponents of the satellite in the u, v, w directions as obtained by super- 
posing the satellite velocities relative to O’ on the motion of that point 
itself). Finally, the Hamiltonian reads 


H = + + Pw?) + pul(vd — *) — p,(ub +76) + V (4) 


and Hamilton's equations for this problem are 


(10) 


The terms V,, V,, and V, represent partial derivatives of the total 
GEm 


potential V = + V with respect to u,v, w. Since 


r=((re + u)? + 0° + 
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“u=—p. + v6 — 7, (8) 
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we have 


le 


( GEm GEm 


( GEm GEm 
— ~ w 
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2 
where we have dropped O (*) ,O (*) etc. 


In the event of no perturbations, that is, u = v= w=0O0 and 
V. = GEm 


7 


, V, = Vy. = 0, Eqs. 7 and 10 become trivial whereas Eqs. 5, 


6 and 8, 9 may be combined to = 7.6? — and 27.6 + 7.6 = 0. 


These we recognize as the equations for conic-section orbits, which are 
satisfied by r, according to definition. Therefore we drop all terms 
from Egs. 3 to 10 which contain 7, but not u, v, w. We may then make 
use of the fact that 


d 
dt 00 
and 
= [GEa(1 — e*) 
where 


a = the semimajor axis, and 
= the eccentricity of the nominal orbit. 


Elimination of p., p., Pw from Eqs. 5 to 10 then yields 


’ — 20’ — u — 2[u + e(u’ — v) sin (6 — 6,) ]/[1 + ecos (6 — 8,) ] 


V./[1 + ecos (@—86,)]* (11) 


v’ +2u’ —v—[—v +4 2e(v’ + u) sin — 6,) + ecos(6 — ] 
+ ecos — 6,)]}* (12) 
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w’’ + [w — 2w’esin (0 — 6,) ]/[1 + ecos (6 — 


where primes denote differentiation with respect to 6, and 6, denotes the 
central angle from node to perigee of the nominal orbit. 

At this point we make the nominal orbit circular, whereby e = 0 
and a = r., the constant radius of the nominal orbit. For the inde- 
pendent variable we shall sometimes write 6 = 6) + f where f = 0 at 
t = to, the starting time of the solution, and 6) locates the moving origin 
O’ at this instant (see Fig. 3). Equations 11 to 13 now reduce to 


— 2v' — 3u (14) 


= = (16) 
2. SOLUTIONS FOR THE RELATIVE MOTION 


Equations 14 to 16 may now be readily solved in the homogeneous 
and the inhomogeneous forms. In the first case we have the solutions 


“= + — + cos f + wo sin f (17) 


2 1 2 
— 3 + f+2 (25, + 3u) 


x sin f (18) 


: wo sin f + wo cos f (19) 


GE 1/2 7 
where 6 = (=) . These solutions are of course dependent on the 


initial conditions for their magnitude and we shall show their applica- 
tion to orbit calculations in the next section. At present we use them 
to generate particular solutions of Eqs. 14 to 16 by a straightforward 
application of the*method of variation of parameters. The results are 
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- |- 2f + 2cosf f cosfaf 


+ cosf f + 2sinf f sin f df 


f f (20) 


—2sinf f sin f df +4.cosf f sin f df 


— 2cosf f fdf| (21) 


w= cos f f f cos (22) 


where all integrals are left indefinite for the time being. 


3. PERTURBATIONS OF THE INITIAL CONDITIONS 


As was mentioned in Section 2, solutions 17 to 19 represent the 
satellite’s motion relative to the circular orbit if at t = tf) some set of 
displacements and velocities uo, vo, Wo, and t%o, do, wo relative to the 
circular orbit exist. These may be errors in the ballistic launch, 
occurring at fo, or observational errors in a sighting which takes place at 
to and from which the orbit is to be predicted. Equations 17 to 19 then 
exhibit the changes of the orbit due to these perturbations in the initial 
conditions and may be used to obtain expressions for the so-called 
orbit sensitivities.‘ 

Since Eqs. 17 to 19 represent deviations of the satellite from the 
circular path, it should be pointed out that the expression for r. + u 
represents indeed an elliptic orbit of small eccentricity. We note that 
the perturbations of initial conditions in a nonrotating coordinate 
system are related to u and v by AV = v» + uo and * = tw) where AV 
and ¢ are the horizontal and vertical velocity errors at cut-off. Solving 
for to and w% from these relations and substituting into Eq. 17 yield 


= + + dr) cosf + (23) 


4 This particular application of moving coordinates to perturbations in the initial condi- 
tions has indeed been made previously by A. D. Wheelon in Lecture 12B of a joint course on 
space technology by the University of California and the Space Technology Laboratories in 
1958. 
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The same result may be obtained by writing the equation of an elliptic 
orbit in terms of its dynamic parameters and expanding for small e. 
We have 


[1 —ecosfcosf, —esinfsinf,] (24) 


r 
where 
cy, = (V+ AV)(r. + Ar), 
cosf, = [(r. + Ar)(V + AV)? — GE]/GEe, 
sin f, = — (r. + Ar)(V + AV)?r/GEe, 


and f, locates perigee from 69. Substituting these expressions into 
Eq. 24, we obtain r = r. + wu as with Eq. 23. 
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Before proceeding to the orbit sensitivities, we note that the com- 
plete set of orbit parameters besides a and e, which have already been 
defined, includes 


@ = the nodal angle between N and y (see Fig. 2), where the node NV 
is the point at which the satellite crosses the equatorial plane from 
south to north and y denotes the first point of Aries as a fixed 
star reference in the equatorial plane, 

= the orbital inclination, 

= the angle from node to perigee in the orbital plane (P denoting 
perigee) where w = 6) + f,, and 
the time of nodal passage or ‘‘epoch.” 
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One now finds it informative to observe the changes produced in these 
parameters by the solution u, v, w. 

The effect of w is to change the location and the orientation of the 
orbital plane. By inspection of Fig. 2 one finds from the geometry 
at the node 


1 1 . 
w(— 90) 4 thy SiN — Wy COS (25) 
sin 2 @ 


and 


Ai = w'(— Oo) = wo COS + Wy sin (26) 


The effects of u, v are of course restricted to changes in orbit shape and 
timing. The orbit attains its apogee at Umax and its perigee at Umin. 
These two points, marking the line of apsides, are seen to lie at 


Uo 


+ 36m)’ 


(27) 


fo.q = 


If 269 + 36u» >.0, a positive wo will put the line of apsides into the 2nd 
and 4th quadrants of f and a negative wp will put it into the 1st and 3rd 
quadrants. Some inspection of the geometry of an ellipse in polar 
coordinates will show that perigee occurs in quadrants 1 and 4 and 
apogee in 2 and 3. With 2%) + 36%) < 0 the situation is reversed. For 
269 + 36u9 = O the line of apsides remains always at right angles to 4, 


where f, = ¥ 3 for w 2 0. Forthe maximum and minimum radial ex- 


cursions from the circle we find 


Umax = + 4 36uo)? + tg? 


min 


from which 


(28) 


2 
Aa = + 


and 
+ 36m)? + wo? (29) 


As should be expected, Aa is independent of wo, that is, insensitive to a 
nonhorizontal ballistic launch. In the solution for v the term v» repre- 
sents an error committed in the original observation of the polar angle 
6) which carries over unchanged into all predictions of passage times and 
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is of no further interest. The secular term in Eq. 18 leads to an in- 
crement of the satellite’s anomalistic period 


Als v(2z) = Vo 26uo). 


It is again worth noting that also AT, is unaffected by tw. The pre- 
diction error for the time of nodal passage is 


= — — %) 
6r.. 


br. [ om 3 (do 26uo)(2nm 60) 
2 + 36u0) sin 2% cos (30) 


where 7 is the number of the revolution for which the epoch is to be 
determined. 

Another use may be made of Eqs. 17 to 19 in designing the corrective 
thrusts required for a satellite rendezvous. Suppose a supply vehicle 
has been brought into the vicinity of some installation that is already 
in orbit and the maneuver of homing in on the target has yet to be 
performed. One may look upon the motion of the cargo ship as a time- 
varying displacement relative to the target of the form of Eqs. 17 to 19. 
Four sets of readings of the displacement vector X(u, v, w) taken by 
the cargo ship relative to the target at four different points of the target 
orbit will furnish the constants in Eqs. 17 to 19. At a suitable instant 


> 
of time ¢;, one may then apply the velocity increment — = and follow 


this up with a second impulse, + z — ¥., at the time ¢, + At where At 


is of negligible magnitude compared with the orbital period. 


4. PERTURBATIONS BY THE EQUATORIAL BULGE 


The well-understood problem of oblateness perturbations shall serve 
as a test case by which to appraise the present method for treating non- 
homogeneous cases. It is customary to make a first approximation to 
the figure of the Earth by introducing the potential 


+ wry? Cir + uy? +0 4+ 


where 7, shall now denote the Earth’s equatorial radius and J is a 
numerical constant. We then find that 
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= (1 — 3 sin? @ sin? z), 


_ 


2 


J 


sin 26 sin?7, and 


In these expressions we have set [(r. + u)? + v? + w?}'? 2r, and 
z~r.sin@sin7z after the differentiations. Insertion of Eq. 32 into 
Eqs. 20 to 22 yields 


Jr2 
c 


[— 1 + sin? + cos 26) ] 


~ [2 — 3sin*i)f + sin 20] 


2 
w= in 27 sin 6]. 


The aperiodic terms in u and v are indeed part of the particular solu- 
tion and reflect the increase of global mass due to the equatorial bulge. 

If we strive for a circular orbit of prescribed size and period, apart 
from the periodic perturbations due to the equatorial bulge, we en- 
counter certain difficulties in realizing this orbit by suitable corrections 
of the ballistic launch conditions. Let us denote the total perturbations 
due to oblateness and launch by @, 3, ®. Then 


Uo + are [— 1 + sin? + 4 cos ] (34) 


c 


i sin 265 (35) 


2 
tio — = sin? 7 sin 265 (36) 


to + [2 + sin? z(— 3 + 4 cos ]. (37) 


The constant term in a may be written 


Ar. = 1 + sin? + 


c 
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and the coefficients of the secular term in @ as 


I Ve 


We may set ip = ti) = 0 right away, thereby imposing conditions on 
vy and ti) which do not interfere with any of the other equations. If we 
are trying to realize a predesigned orbit, we should certainly wish to 
have the mean values of r and 6 agree with the specifications and there- 
fore demand that the right-hand sides of Eqs. 38 and 39 vanish. Un- 
fortunately, these equations are not algebraically independent in » and 
uy. We may adopt one or the other and then take Eq. 34 or 37 as an 
additional condition. 

If the emphasis is on an orbit of the correct mean radius, we may 
take Ar, = 0 and ad = 0 where the latter condition calls for cut-off at 
nominal radius. Solving Eqs. 34 and 38 for u» and #» and substituting 
the results into Eqs. 37 and 39 we obtain 


iy = 1 + sinti( + 26,)] = AV (40) 


c 


2 


Aé 


If on the other hand we take Aé = 0 and 


we obtain 


(42) 


(43) 


Equations 40 and 42 vanish only for certain combinations of 7 and 6 
while Ar. = Aé = 0 for 7 = sin— v2. 
Since the sin @ term of w in Eq. 33 represents a perturbation on 14, 


we may use Eq. 19 to eliminate it. The resulting corrections for the 
ballistic launch are 


and 
6Jr2 


sin 27 cos 4». 
2r. 


and 
AV = to + tof = 0, 
Uy = [3% + sin? z(— 1 + 3 cos 26») ] 
and 
Ire 
Ar. = — — } sin*¢]. 
re 
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St, 
4.23 


As a numerical example let us take r, = 5000 statute miles, 7 
and 6) = 50°. Then Eqs. 40 to 44 yield AV = — 4.37 ft/sec, A6 
10-8 rad/sec, i) = 245 ft, Ar. = 858 ft, and wy) = 7.65 ft/sec. 
The periodic terms of u and v in Eq. 33 have the double angle 26 
as their argument because of the biaxial symmetry of the force field 
produced by the equatorial bulge in the orbital plane. Since the attrac- 
tion from the bulge acts in opposite directions on the northern and the 
southern halves of the orbit, the simple argument 6 is to be expected 
in the periodic terms of w. The term fcos@ in w indicates a pre- 
cession of the orbital plane and the nodal regression is found from 


Aw = “ as 


Aw (GE)? GE 
= Aw ae = — Jr2 cost 
sin 2ar,3/? rf 


which agrees with the results obtained by other methods and amounts 
to 5.80 10-7? rad/sec or 79 = 125 days for the nodal period of a 5000- 
mile orbit with 7 = 51°. Finally one notes that the complete expres- 
sions for u, v, w with the initial conditions adjusted as above furnish a 
convenient means for predicting satellite positions throughout the orbit. 
The maximum amplitudes of the periodic perturbations impressed on a 
nominal orbit of 5000-mile radius are approximately 7000 ft. 


5. PERTURBATIONS BY EXTRA-TERRESTRIAL MASSES 


For a study of gravitational perturbations by the moon or sun we 
locate the perturbing body P in the x, y plane where the latter does not 
in general have to contain the equator. The central angle 6 = 6) + Af 
locates P in the x, y plane in a fashion analogous to the way in which O’ 
was located in the orbital plane (Fig. 3). 4 is measured from the posi- 
tive x-axis. 6 is its initial value at ¢ = f = 0 and X is the ratio between 
the angular velocity of P and the satellite. We note that the present 
problem will be rheonomic not only because of the time-dependent u, v, w 
coordinate system but also due to a time-dependent perturbing potential 
resulting from the motion of P. \ will be < 1 in all cases of practical 
interest. 

If one writes the perturbing potential as 

1 XX» + + 


where Xp, 2p are the coordinates of P, R,? = + y,? + and 
Rom? = (Xp — x)? + (vp — y)? + (2, — 2)? then one finds for the per- 
turbative acceleration of the satellite in the x-direction 


: 
Fa 
= 
y 
m ax R,' 
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0X, Y,2) 


x 
Fic. 3. 


which is the sum of the direct and indirect perturbations on the satellite. 
Analogous expressions follow for 7 and 7. We use approximations of 


the type 
mGP | 3x, 


(xxp + + — 


and take x y ~r.sin@cosi, ~r,sin@sin7z in these ex- 
pressions to obtain 


— [4 cos?i — sin?¢ + 2 sin?z cos 20 + 2 sin? 7 cos 26 


+ (1 — cos 7)? cos 2(6 + 9) + (1 + cos 7)? cos? (6 + 6) ] 

&P,* 

+ (1 + cos 7)? sin 2(6 — 6) — (1 — cos 1)? sin 2(6 + 6) ] 


[— 2 sin? sin 29 (45) 


3r.GPm 


= [sin 27 sin 6 + (1 — cos sini sin (26 + 8) 


+ (1 + cos?) sini sin (26 — 6) }. 


t 
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Executing formulas 20 to 22 one finds 


where g = 


the case \ = 0 separately, formulas 20 to 22 yield 
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“u= - cos? — sin? 2) 

2 

+ cos 26 + sin? cos 26 
(A + 2)(1 — cos i)? 

(A + 1)(2d + 1)(2d + 3) 


(A — 2)(1 + cos 7)? 
(A — 1)(2X — 1)(2d — 3) 


cos 2(6 + @) 


v= 402 cos? i — sin? 7) f 


sin? i 
d (44? — 1) 
+ 12d + 11) (1 — cos i)? 


(442 — 12 + 11)(1 + cos 7)? 


HA = 3) 


— 4} sin?¢ sin 26 — sin 26 


9 E sin 27 sin 6 — f sin 27 cos 6 


(1 — cos?) sinz 


hast) sin (26 + @) 
(1 + cos 7) sini 
— 1) 
3Pr.4 
8ER,* 
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cos 2(6 — | (46) 


sin 2(6 — | (47) 


sin (26 — | (48) 


This solution is invalid for } = 0, +1, +3, + 3%. Treating 


u = gL— 3(2 cos*z — sin?) + $sin®icos 26 — 2 sin*i cos 26, 
+ 3(1 — cos7)? cos 2(6 + 60) + + cos7)?cos2(@ — (49) 
v = gL 4$(2 — + 3 sin’? cos 26)) f — 4. sin? 7 sin 26 


— (1 — cos7z)? sin 2(6 + 60) — + cos i)? sin 2(@ — 6) (50) 


- 

| 
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w= gf (1 + cos 7) sini cos (6 — 
— (1 — cos?) sinz cos (6 + 26)) — sin 27 cos 6 ]f 
+ 3} sin 27 sin 6 + 3(1 — cos7) sinz sin (6 + 26)) 
— }(1 + cos 4) sin isin (@ — (51) 


For ease of physical interpretatioa, we begin with an examination of 
this case. 


If we set x, = 2, = 0, y, = Ry, bo , t= 0, 6 = 0, and f = 0 


then we have coincidence between the orbital planes of the satellite and 

the perturbing body; in addition, the latter has been fixed on the posi- 

tive y-axis and the satellite is assumed to start its motion on the positive 

x-axis. With these conditions we find 

37.GPm 

3r.GPm 
sin 26 


[3 — cos 26] 


= — $ cos 26] 
+ 44 sin 26] (49’) 
= (), 


The first terms in u and v are due to the first term in V, of Eq. 45’. 
The secular term in v, which is to be avoided in the case of a 24-hr 
satellite, may be eliminated by use of the complementary solutions in a 
fashion similar to the calculations of Section 4. 

The second terms in « and v are consequences of the so-called “‘tidal 
effect”’ of P on the orbit. One notes from u that the orbit elongates in 
the x-direction and flattens in the y- direction and 9 reflects the angular 
anomalies which the satellite must undergo in order to satisfy conditions 
of angular momentum and energy while its distance from the central 
mass is continually changing. In fact, these considerations of angular 
momentum and energy permit an alternate way of analyzing the present, 
simplified case which will be reproduced here for illustrative reasons. 

Figure 4 shows a view of the x, y plane and the forces, due to P, 
acting on the satellite in its orbit. It can be shown that 


cos 0 


V. = 0 
and 
v 
w 
Up 
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and 
2r.4Pm sin 6 cos 4 
= 
R,? 
and thus, with z = 0, 


3r.GPm sin 6 cos 6 


3 


[2 sin? @ — cos? 6]. 
Pp 


We recall that 
dK 
“dt FLU, + F,U, 


dA 
Fi(r. + u) — Fy 


= total energy of the satellite in the geocentric frame, 
angular momentum of the satellite about the Earth, 


U, = u — v6 
U, 6(r. + u) + 3, 


and 
a | 
= 
/ 
/ \ 
/ 
\ 
Fy 
Fic. 4. 
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where 
and 
i 
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the instantaneous geocentric radial and tangential velocities, respec- 
tively. From 


GEm 
+ u)? + 


m 


A = + u) — (u’ — v)v] 


we obtain after simplification 


E 
[r.(2u’ +0) + u(2u’ +0") + 0(20’ — + — 3u) + 


= [3 sin + u + v’) + (2 sin? 6 — cos? 6)(u’ — v) 


and 


Lre(2u’ + + u(2u’ + + (20' — 
_ Pr. 
R,' 
To O(r.u) inside the parentheses on the left and to O(r.?) on the right, 
both equations yield (15) and to O(u*) on the left and O(r,«) on the 
right one may combine them to obtain (14). 
In a case where 7 = 0 and } # 0, + 3, + 3—for example, \ ~ 35 
which is the ratio of angular speeds between the moon and a 24-hr 
satellite, we get from Eq. 46 


i= 


[3 sin 6 cos + u) — (2 sin? — cos? 6)v }. 


cos 2(1 — (40’) 


(2n + 
2(1 — A) 
this means that the orbit’s minor axis, generated by the tidal effect, 
rotates in phase with the perturbing body P. The periodic term in 7 is 


Umin OCCUrS at 6 = ,n = 0,1, 2,---. As was to be expected, 


(442 — 12d 4+ 11) 


which produces maximum circumferential displacements of the satellite 
at points which are 4 


h = gy we find Yuax ~ 3000 ft for a 24-hr satellite. If the maximum 
angular displacement of such a relay must stay within + 10° for trans- 


out of phase with respect to the orbit axes. For 


t 
and 
= 
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mission purposes, the periodic term which was just computed is seen 
to be negligible. If, however, one assumes a typical guidance error of 
io = AV = 1 ft/sec during the launch, then the time between correc- 
tive tangential impulses for such an installation will be of the order of 
10r.6 
57.3967 
Let us next consider a case where z # 0, but \ = O and all other 
previous conditions are maintained. The tidal distortion of the orbit 
is now somewhat more involved as can be seen from Eqs. 46 and 47 
and may be understood by inspection of the geometry. However, the 
perturbing body now also exerts a tidal moment on the orbit which 
should give rise to precession. Indeed we find from Eq. 51 that 


= 93 days, where we have used Eq. 18. 


w = gsin 2i[sin 6 — 26 cos 6] (51’) 


where we shall limit our attention to the periodic term of linearly in- 

creasing amplitude in 6. One observes that wmax in a given cycle due 

to this term occurs at the node and that its rate of change as a function 
AWmax 3 sin 2¢Pr,4 


of @ between epochs is  4ER,* where is the increment 


iN Wmax per revolution. 


A 


TIDAL FORCES ON THE 
ORBIT (REVERSE SENSE FOR 
RADIATION MOMENT) 


Fic. 5. 


On the other hand, arguing again from fundamentals, we have with 


M ‘ 
the help of Fig. 5: ¢ = oe where M is the average tidal moment 


exerted on the orbit (positive in the right-hand sense about the x-axis), ¢ 
the average angular rate of orbital precession about the y-axis resulting 


ts 
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| | 
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from this moment (positive in the right-hand sense) and A = m(GEr,)!/? 
is the average angular momentum of the satellite. From the geometry 
at the node we have ®Wmx = — ¢f-cost. Now, 

1 


Jo 


(F,r. sin 6sint — F,r, sin 6 cos 1)d6 


M 
where 
rGPm sin 6 sin 7 

R,? 


resulting in 
3 sin 
4R,, 3 


Combining all of the necessary relations, we have 


M= 


r. COS 4 1 3 sin 2ir2GPm 3 sin 2ir,4P 


m(GEr,)'/? cos 4R,' 4ER,? 
The most relevant effect of this precession, as far as orbit acquisition 


is concerned, lies in the nodal regression 


cos 4P(GEr,*)!? 
2m sin ir, 2ER,' 


Since this expression assumes a stationary perturbing body P it will 
be more realistic to consider the results for \ # 0. From Eq. 48 
we find 
= 3 cos tP(GEr,’)'/? 

4ER,! 
There are two remarkable features about this result: (1) it is inde- 
pendent of @ and (2) it is one-half the precession rate with stationary P. 
As } — 0 one may show that the last two terms in Eq. 48 degenerate to 
a cyclic term of increasing period which in the limit contributes another 
secular term to account for = 

Let us assume that the lunar orbit and the ecliptic coincide with the 
equatorial plane® and let us suppose that we have a 24-hr satellite with 
an orbit inclination of 7 = 20°. Then the nodal regressions due to the 
moon and the sun are Qy = — 0.82 10-* rad/sec and 2, = — 0.30 10-9 
rad/sec which add to 2, = — 0.260 10-§ rad/sec from the equatorial 
bulge. The decrement in nodal period due to these sources can be 
written as 

27 (Q M 2,) 
+ Qa + 
5 Without this assumption, the orbit inclination with respect to the x, y plane (contain- 
ing P) would be variable as a result of the equatorial precession. 


AT. = — 8418 days. 


| 
ome 
= 
4 
‘ : 
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(T, due to oblateness only is 27,956 days.) The decrement in geo- 
period (time between successive ascents over the equator)® may be 
written as 


— 1.24 sec 


for the moon and sun combined. If a secular perturbation of 7 due to P 


— dt 
were to exist, it would follow from 7 Aw’(@ = 0) but this is seen from 


Eq. 51’ to vanish, as should be the case for small rotations of the orbital 
plane about the y-axis. 


Fic. 6. 


6. PERTURBATIONS BY SOLAR RADIATION PRESSURE 
The total force exerted by the sun’s radiation on the satellite may be 


written as pe where 6 takes into account the shape, size, reflectivity, 
pm 

and mass of the satellite and R,,, is the distance from the sun to the 

satellite. Since the pressure force obeys the inverse square law, we may 


write V = B so that F; = — is the pressure component in a given 


Kon 0g 


direction. This represents the radiation as coming from a point source. 


6 Also known as draconic period. 
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Its divergence is unrealistic in the case of solar radiation in terrestrial 
space and will be eliminated at a later point. In general the illumina- 
tion of the satellite orbit will be nonuniform due to the Earth’s shadow, 
albedo, and atmospheric absorption. We represent this nonuniformity 
as follows. Assume that the radiation source is fixed in space—then 
the orthogonal projection of the line EP on the orbital plane marks the 
constant angle 6) from the node (Fig. 6) and we introduce 6 = f + 4 
— 6) = f + 6. Now the distribution of radiation intensity around the 
orbit will be approximated by the first two terms of a Fourier cosine 


series so that V = (1 + bcos 6). With the usual simplifications 
we get 


Bmr,. 


= R <[- 1 — bcos 6+ £cos (gcosi + Zsin7) sin 6 
Pp 


+ + Zsin2)[sin (6 + 6) + sin — 60) 
+ (6 + 6) + cos — 0) J], 


Bmr,. 


— 3bz[sin (6 + 6) + sin (@) — 
+ 4b(g cos i + zsin i)[cos (8 + 4) + cos (6) — 60) J], 


[ Esin 6 + (7 cost + Zsin 7) cos 6 


= (Zcost — gsinz)(1 + 6), 
Pp 


where = —, 1 


Substitution of these expressions into 


Eqs. 20 to 22 yields 


“= — cos — 60) — 4(gcosi + Zsin2)b sin (69 — 
+ f[bz sin — 60) — b(gcosi + Zsin 7) cos (6) — 
+ 3bsin 6 — + + Zsin72) cos 6] 
+ {bcos 6 — f#cos 6 — + Zsin7) sin 


+ cos (6 + 6) + 3b(gcosi + i) sin (6 + 6)] (53) 


| 
rye 
: 
hi 
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v = 2bz sin (6 — 60) + b(g + Z sin 4) cos (6) — 
+ f[#(b cos (00 — 60) — 3.cos 6) + (gcosi + Zsin 7) 
X (b sin (60 — 4) — 3 sin @) — 2 + bcos 4] 
+ (9 cos + Zsin 7) cos — 60) — — ] 
— dbsin 6 + 6 — cosi + Zsin 7%) cos 


— Hz sin (0 + 4) + cosi + Zsin i) cos (6 + (54) 


w = 2(gsini — Zcos7)(1 + }dcos 6 + sin 6) (55) 


GER,* 

We examine these results again for some special cases. If x, = y, 
= 0,1 = 0, b = 0, zg, = R,, the orbit lies in the x, y plane and is fully 
illuminated from the positive z-axis. Then u = 3, v = — 2f%, and 
w = 2g. The meaning of w is obvious. It represents a slight transla- 
tion of the orbit out of the x, y plane due to the uniform radiation pres- 
sure acting on it. From Fig. 7 it is clear that for m to be in equilibrium 


where 7 = 


4 397 
\ 

| 
\ 
| 
| ia 
| 
Yon le 
| Fic. 7. 
| 


398 F. T. GEYLING (J. F. 1. 


one requires — = GEm rom which w= — 
time-independent “ component represents an enlargement of the orbit 
due to the radial pressure arising from the divergence of radiation from 
P and the linearly decreasing v is then necessary to conserve angular 


momentum. From the last condition 


(6 + Ad)(r. + u)? 


(6 + Aé)? (1-44+0(4)). 


From radial equilibrium 


= 27g. The 


and thus 


Bmr, _ GEm 
(re + 
GEm 


m(6 + Aé)?(r. + u) + 


and with the above expression for (6 + A6)?, we have after omission of 


u* 
O { — } the desired result for vu. From the momentum equation one 


finds Aé = 2—é6and thus v = — 2fz. 

With x, = z, = 0, y, = R,, 1 = 0, 6 = 0 we now turn to the hypo- 
thetical case of an orbit that is uniformly illuminated by a source in its 
own plane. ‘This example is a worthwhile stepping stone to the under- 
standing of the same case with nonuniform illumination. From Egs. 53 


and 54 
2039 cos — + 1] (53’) 


— o[39 sin + 2] — 4A9 cos (54’) 


dA 
In order to interpret this result by means of ——- and —— as explained 


dt dt 


in Section 5, we observe that now F, = — — sin @and F, = 


R;? 
as seen from Fig. 4, if one neglects the small effect of the divergence of 
radiation from P. Going through this elementary approach, one finds 


u = cos@ — 
and 
— sin + 4 cos 6] 


Gs 
u u? 
ag 
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which differs from the exact expressions by the small terms 7 and — 262 
Bmr, 
substitution into Eq. 14. This discrepancy is, of course, due to the 
assumption of perfectly collimated light, which is true for solar radiation. 

We observe that the amplitudes of u and v in this solution grow 
without bound as @ increases. One would indeed expect this result 
from the linearized equations for small oscillations about the circular 
path if the exciting forces F, and F; occur at resonance with the orbital 
frequency—which is the case. In terms of orbit shape this linearized 
model predicts ever-increasing eccentricity. Solving for the line of 


du 


apsides, we have from = 0 that 


and reflects the omission of the term — in V,, as may be seen from 


(56) 


6= — §ctné. 


This results in the perigee locations 


w, = 3.08 + (3.08) 3.125 + + 5.425) 


for successive passes, where ,_ is a very small correction. (A corre- 


e 


sponding sequence of roots may be written for w,, the exact apogee 
location.) Now 


Umin = + w,) cos wy — 7 sin w, |, 
and (57) 


Umax & a(x (2n +1)+ COS + F SiN w, |, 
n= 0,1,2---. 


If we let the average value of the semimajor axis, a, throughout one 
revolution be defined as 


a=f, + (Umax + Umin) = COS 


then it becomes apparent that 2a of the orbit remains unchanged. For 
the average eccentricity, however, we get 
1 


é= (Umax — Umin) 6(x(2n +3)+ cos w, + 7 sin w»| 


fiers: 

and 

= 

dé 329 

~~ 

— = = — ~~ COS w». 

dé 2r. 
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For a 100-ft mylar balloon traveling in a 5000-mile orbit we find 


p= = 1.65 10-/sec. 


In a period of one month this would result in an eccentricity of 0.428 
which is indeed prohibitive. Long before this condition is reached, 
other physical phenomena interact with the radiation pressure to shape 
the orbit. A more complete discussion of these perturbations will be 
given in Section 7. 

If the orbit illumination is nonuniform, 6 # 0, then we have in 
place of Eq. 56 


6= — ictné + $b cos (56’) 
With } = 1, this yields 


v(3.10 3.126 


c 


In place of Eq. 57 we obtain 


Umin = BL + w,) cosw, — w, — b(} — cos | 


and 


Umax & — 3x(2n + 1) cosw, + sin w, — b(4 + cos w,) 
n=0,1,2---,] 


where the absence of a secular term in 6 happens to be due only to the 
fact that 6 = — 4 From this one can find estimates of d@and é similar 
to the ones for 6 = 0 and subject to the same limitations. 

We finally take 7 # 0 and maintain x, = 2, = 0, yp = Rp, 00 = 9, 
and b #0. Then 


w = + 3b sin 6 — 63) cos 6). 


The secular term in this result again represents a precession and is due 
to the radiation moment M exerted on the orbit due to the term b cos 6 
of the intensity distribution. Using the arguments of Section 5, it can 
be shown that 


bBm sin tr, 


M = 
2R;* 


4 
ext 
(57’) 
te 
: 
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From this or from the expression for w one finds 


b sin 2 (58) 


If we take a = 5.66 10-‘ ft/sec? for the 100-ft mylar balloon, then 
Pp 

1.228 rad/sec and = 5850 days for ; 

= oRaNGE = 1, rad/sec and T, = 5850 days for a 


5000-mile orbit due to radiation pressure alone. The increment in 
the geo-period is AT¢ = 0.095 sec if 7 = 20°. One also observes that 
“ = 0, as in the case of Section 5. 

As a last numerical result, we estimate Az as produced on a 5000-mile 
polar orbit during one-eighth of a year when the sun travels in the x, y 
plane from the positive x-axis halfway to the positive y-axis (Fig. 6 
with 7 = 90°, 6) = 0). We note from Eq. 58 that the precession due to 

AW mind 
about the 
COS 
line EP. Reinterpreting this result for the present situation we find 
dit _ Br.’ sin a6 
dt 2GER;? 
Thus, 


the radiation moment consists of a rotation ¢ = 


where a is the angle from the positive x-axis to the sun. 


where o, is the angular rate of the sun in the geocentric frame. During 
the 2nd and 3rd octant of the sun’s orbit the satellite is assumed uni- 
formly illuminated and no precession takes place. In the 4th octant 
an equal and opposite effect to that of the 1st octant occurs. Thus the 
net effect of precession under radiation pressure vanishes in each half 
year. 


7. COMBINED RESPONSE OF THE ORBIT TO RADIATION PRESSURE 
AND THE EQUATORIAL BULGE 


In the preceding section we noted that the linearized equations of 
motion for a satellite on a nominal circular orbit give rise to unbounded 
u and v in response to radiation pressure within the orbital plane. 
At first sight one might expect the nonlinear equations of motion for 
this problem to impose a restriction on the magnitude of the response in 
the perturbations—in particular, one would hope for a bound on the 
eccentricity which keeps the perigee of the orbit at a sufficiently great 
altitude to avoid orbital collapse from atmospheric drag. In order to 
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perform the nonlinear analysis, we now include terms to O { — } inthe 
expressions for the attraction by the central mass. Thus, 


Ou r 


27 


GE G 7D 3 6 
ml =; uv 


r r3 


3 ut 
tee) +0(%)| 
and 


0 {GEm 3 6 


If we substitute these expressions into Eqs. 5 to 7 together with 


+ uw?) +O 


m 
cos 6, and V,, = 0 (since we have = 0, 


= 0, and y, = R,), then we may again eliminate p,, p., Pw from 
to 10. We obtain 


= 
Eqs. 
E, = — 20’ — 3u + bu? + b.(v? + w?) 

+ + b4(uv? + + Visine =0 (60) 
E, = + 2u’ — + + b.(v* + vw?) + V; cos = 0 


where 


GER? 

and £, and EF; shall be used as shorthand notation for the differential 
operators on « and v._ By inspection of the linearized solution of Sec- 
tion 6, we may guess at the form u = u, sin 6 + u.cos @andv = 9; sin 6 
+ v.cos 6. We now use this form to develop an approximate solution 
of Eq. 60 by the Ritz-Galerkin method. The Ritz conditions to be 


Age ov 
: 
w=0 
3 3 27 6 
b= or. = or?’ by = 2’ 
j6 
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satisfied are 


f E,cinede-= 0, f 
0 0 


** Fy sin = 0, f 
0 0 


Performing these quadratures, we obtain the four nonlinear algebraic 
equations in Us, and v2: 
4u, — 2v. — + 

— + + = V; 

— + = 0 

+ — + bov,*) 

+ + 3b.0,02') = 0 
— 2u, + v2 (b.(2u + + 3b.0,*02) 
— + by.*) = Vi. J 


In the normal course of events one would now solve these nonlinear 
algebraic equations by a suitable numerical iteration scheme—and 
presumably a similar set of equations would be written and solved for 
b #0. Unfortunately, however, the nonlinearities governed by the 
coefficients 6, to b, in this problem turn out to be so feeble that the 
dynamic response of the satellite is very nearly described by the results 
of Section 6 if no additional phenomena enter the picture. 

It turns out, however, that precession of the apsides due to the 
Earth’s oblateness has a decisive influence on orbit distortion in this 
case. The reader will be aware of this shift in perigee from the litera- 
ture. It is a nonuniform motion for inclined orbits but for equatorial 
orbits it simplifies to a steady rotation of the line of apsides, at the rate 


= 2f (62) 


Wp 


where positive w, is taken in the sense of the anomalistic motion ). 
This precession occurs if an eccentric orbit, once it has been established, 
is left under the influence of the Earth’s oblateness only. 

Now let us recall the linearized effects of radiation pressure in a 
simplified form before we combine them with the perturbation from the 
equatorial bulge. We have learned that under the action of radiation 
pressure alone the secular term in u distorts the orbit according to 
i = Cé cos 6 where perigee lies at 6 = 7 and apogee at 6 +0, C = $39, 


: 
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(and the symbol @ is used in distinction from u which will be reserved 
for the composite effect of oblateness and radiation pressure). This we 
have shown to be very nearly true after a small number of revolutions 
with the help of Eqs. 53’ and 56. Now it must also be pointed out 
that the orbit distortion as given by @ does not exhibit the exact form 
of a cosine function for any one complete revolution. This is due to the 
time-dependent amplitude term in #. The cosine shape describes an 
osculating locus only which the satellite would follow if at any one 
instant the perturbing radiation were eliminated.?. If however we 
assume that C < 1, then the individual pass of the satellite very closely 
approximates a cosinusoidal departure which is superimposed on the 
nominal circular shape but whose amplitude changes appreciably only 
over a large number of revolutions. This assumption is a convenient 
one and a reasonable one to make in the present case. 

Now we proceed to combine what we have learned about the indi- 
vidual actions of the equatorial bulge and the radiation pressure to 
obtain the following picture. During each revolution of the satellite 
a distortion z, = Cf7 4 cosf is impressed on the orbit by the radiation 
pressure; but in the same time all such terms that were previously 
generated have shifted along the orbit by the phase angle Ay = #,7 4, 
where 7°, is the anomalistic period. The composite perturbation of the 
orbit is then 

Ni 
F(@, N:) = ¥& CfT 4 cos (@ — ndy) (63) 

n=0 
up to N, revolutions of the satellite. In writing this expression we 
have not only applied linear superposition to the perturbations from 
oblateness and radiation pressure but also to the effects of individual 


revolutions. Substituting 7, ad and ¢ = = and letting Ay — 0, 
Wp W» 


N, — ~ such that N,Ay = 7: we have 
Ni 


F(0, yi) = lim >} cAy cos (f — ndy) 


Ay0 n=0 


cos (f — y)dy = c[{sin @ — sin (6 — ] 


= 2¢ sin cos (0 ). (64) 


“ 


’ Another way of visualizing the osculating orbit may be obtained by the following ‘‘Gedan- 
ken-experiment": Imagine that the satellite is dissolved into a string of particles continuously 
distributed around the orbit (for example, an artificial ionosphere). Then, barring all com- 
plementary solutions by a suitable choice of initial conditions for each particle, this ring of 
distributed mass will indeed exhibit a cosine distortion at every instant whose amplitude how- 
ever is time-dependent according to @, 


= 
Bes 
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We see now that the resultant orbital distortion is again a cosine 
function in 6 which, however, proceeds around the orbit at an angular 
rate which is half the precession due to oblateness. ‘Thus, the apsidal 
motion of the resultant orbit is f, = }a,. Since the amplitude of 


the perturbation varies periodically in Had we note in particular that 


2 
Fin = — 2c sin Thus, perigee first originates at f, = 7, then shifts 
3r 
tof, = > where it attains a maximum amplitude of —2c, proceeds to 


f, = 2m where it vanishes, and immediately reappears at f, = 7 since 
what used to be apogee now turns into perigee with the change of sign 


in sin of Eq. 64.° 

One notes that no changes need be made in the above analysis if 
the line EP were to rotate with respect to an inertial reference frame at 
the rate o,. We merely change the notation from @ to 6, and replace 
@, by a, = w» — o,, to indicate that these quantities are measured in 
the rotating frame. 

If we take the maximum perigee depression as a measure of orbit 
feasibility, we have 


For a 5000-mile orbit this amounts to 115 mi. and for a 26,250-mile orbit 
to 11,400 mi. We also have for the 5000-mile orbit: #,, = 1.6 10-® 
rad/sec; duration of one sweep of perigee = 45 days; and range of one 
sweep of perigee = 3.91 rad (in an inertial frame). For the 24-hr orbit: 
®pr = — 1.93 10-7 rad/sec; duration = 377 days; and range = 3.33 
radians. 


8. SUMMARY AND CONCLUSIONS 


Reviewing the material in Sections 3 to 7, we note that Hamilton’s 
canonical equations lend themselves remarkably well to the solution of 
satellite perturbation problems in terms of position coordinates versus 
time. Orbit sensitivities to errors of injection and observation as well 
as the perturbations from various geophysical effects may be calcu- 
lated. We recall in particular the following numerical results from the 
analyses of circular orbits. 

If we wish to establish an orbit whose shape after discounting 
periodic perturbations due to the equatorial bulge is a circle of 5000-mile 
radius at an inclination of ¢ = 51°, with the cut-off point at 6) = 50° 


8 It is also seen that Eq. 64 gives the shape of an artificial ionosphere as a function of yi, 
that is, time. This does not include “‘tidal effects.” 


: 
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from the node and at a cut-off radius of 5000 miles, then the necessary 
change of injection velocity from the circular value is AV = — 4.37 
ft/sec and the correction in azimuth Ay = — we ~ — 0.33 10-* rad. 
In spite of these corrections the angular motion will show the dis- 
crepancy Aé = 4.23 10-8 rad/sec. If we insist on a mean velocity and 
launch velocity equal to the circular value, we need to cut off at an 
increase of altitude a = 245 ft and the resultant mean radius will be 
increased by Ar, = 858 ft. The maximum amplitude of the periodic 
perturbations which are impressed on the circle amounts to about 
7000 ft. 

The periodic perturbations on a 24-hr satellite in the equatorial 
plane by the moon’s gravitation will amount to about 3000 ft in the 
circumferential direction which, indeed, is negligible for practical pur- 
poses. However, a guidance tolerance of +1 ft/sec in the injection 
velocity or subsequent corrections of circumferential velocity to that 
accuracy will ensure an operating period of only 93 days before the 
satellite drifts out of an allowable geocentric sight-angle of +10°. 
If a 24-hr satellite were put into orbit at an inclination of 20°, the regres- 
sion of the nodes due to the moon and sun combined would amount 
to 2 = — 1.12 10-* rad/sec. This reduces the nodal period of 27,956 
days, which is due to oblateness alone, by 8418 days. The decrement of 
the geo-period amounts to — 1.24 sec. 

Solar radiation pressure acting alone on a 100-ft mylar balloon of 
0.0005-in. thickness produces the following perturbations. If the or- 
thogonal projection of the line ES on the orbit plane marks the angle 
6, from the node in the direction of satellite motion then a perigee 
will be impressed on the initially circular orbit at approximately 


69 + 2 This location of the apsides is quite insensitive to the inten- 


sity distribution of the incident radiation around the orbit. Further- 
more, no appreciable change occurs in the semimajor axis of the orbit. 
The rate at which eccentricity is being accumulated, however, is of the 
order 1.65 10-7/sec for a 5000-mile orbit. According to this analysis 
the satellite would scarcely have a lifetime of a month. If, however, 
the precession of apsides due to the equatorial bulge is included, then 
the following perturbations result. A 5000-mile orbit in the equatorial 
plane will develop a periodic perigee with a maximum depression of 
Umin = — 115 mi. and a precessional rate of f, = 10-* rad/sec in an 
inertial frame. The duration of one sweep in perigee between zero 
amplitudes is 7, = 45 days and the angular range through which the 
perigee travels in each sweep relative to an inertial frame is A@ = 3.913 
rad. The corresponding quantities for a 24-hr satellite are: Umin 
= — 11,400 mi.; f, = 1.03 10-7 rad/sec; T, = 377 days; and A@ = 3.33 
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rad. If amylar balloon were put into a 5000-mile polar orbit, then the 
moment exerted by radiation pressure on this orbit may perturb the 
inclination by as much as 6.2 10-* rad in one and a half months—the 
net effect over each half-year, however, reduces to zero. If the same 
orbit is inclined at 7 = 20°, we obtain a nodal precession rate of the 
order 1.228 10-8 rad/sec which turns out to be 2, in magnitude of the 
regression due to the equatorial bulge for an orbit of this size. The 
increment in the geo-period amounts to 0.095 sec. 

Some order-of-magnitude checks of the above results have been 
obtained by a digital simulation. Work in this direction continues and 
will be reported at a later time. 
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A FURTHER CREEP ANALYSIS FOR ROTATING SOLID DISKS 
OF VARIABLE THICKNESS 


BY 


B. M. MA! 


ABSTRACT 


A further creep analysis for stress distributions in rotating solid disks of variable 
thickness and at uniform temperature is presented herein, based on the Tresca cri- 
terion and its associated flow rule. The advantage of this analysis lies in the fact 
that the complicated problem is to be solved in great simplicity. 

Disks of gas turbines and jet engines are usually operated at relatively high 
stresses and high temperatures. Therefore, in this paper the stress analysis for the 
variable-thickness and uniform-temperature disk at steady-state conditions will be 
analyzed and discussed by using the exponential function stress-creep-rate relation. 

Illustrative examples for computing the stress distributions of the variable- 
thickness and uniform-temperature disk will be presented on the basis of theoretical 
analysis. For the purpose of comparison the stress distributions for the constant- 
thickness and uniform-temperature disk will also be given. It will be found that the 
stress distributions over certain central portions of the disk are quite different between 
the variable- and constant-thickness disks. 


NOMENCLATURE 


The following nomenclature is used in the paper: 


a = constant 
c, m = constants 


h = f(r) = cr™ = cr.™x™ = thickness of disk 
k, K = constants 
n = material creep number 
r = radius to any point of disk 
r; = inner radius of disk 
ro = outer radius of disk 
r. = radius to the point x, of disk 
t = time 
T = absolute temperature 
it = rate of radial displacement of disk at radius r 
x=r/r, 
= ri/To 
Xe = fe/Te 
o« = equivalent or effective stress 
= pw*r,* = equivalently inertial stress in rotating thin ring of radius r. 
= material creep modulus, having dimensions of stress 
72, 7; = principal stresses along axes 1, 2 and 3, respectively 
o,, 0, = tangential, radial and axial stresses of disk, respectively 
o,, = radial stress at r, of disk 
o,, = radial stress at x, 
= tangential stress at x, 
creep rate or effective creep rate 
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1, €2, €; = principal creep rates along axes 1, 2 and 3, respectively 

€:, €r, €: = tangential, radial and axial creep rates of disk, respectively 
€1, €2, €3 = principal creep strains along axes 1, 2 and 3, respectively 

p = y/g = density of disk material 

y = specific weight of disk material 

g = gravitational acceleration 

w = angular speed of disk 


1. INTRODUCTION 


In a previous paper (1),? some formulas were derived for calculating 
creep deformations and stress distributions in rotating solid disks of 
variable thickness and uniform temperature. These formulas were 
based on the maximum shear theory associated with the Mises flow 
rule, and also under the restriction that the tangential stresses are in- 
variably greater than the radial stresses in the disk, except at the center 
of the disk. The results obtained appear to be relatively simple, reason- 
able and in excellent agreement with the creep test data (2) when com- 
pared with the results found using the distortion energy theory. 

This paper, however, is based on the theory of the Tresca Criterion 
and its associated flow rule. The restriction on the stresses a; > ¢, 
over the whole disk is also removed. In general, the stress distribu- 
tions in the disk are divided into the two regions below: 


l1.o:>0,>0 for Se 
2. =o, > 0 for 


and the stress equations are obtained in the closed form without re- 
sorting to the successive approximation method. 

Stress distributions in region 2 of the disk can be analyzed by recog- 
nizing the fact that radial stresses are frequently equal to tangential 
stresses near the center and over a portion of the disk (3). This is true 
because the centrifugal loading or inertial forces of turbine blades are 
transmitted through the rim of the disk to the disk proper, that is, the 
radial stress o,,, acting on the outer radius of the disk. The greater 
the value of ¢,, is, the more the stress condition, o, = a,, of region 2 
over a portion of the disk spreads toward 7,. 

Experimental data also show that the yield and bursting speeds of 
hollow disks are materially reduced by the presence of a bore, because 
the presence of this bore is equivalent to introducing a stress concentra- 
tion factor at this point of the disks. The reduction of these yield and 
bursting speeds is generally in proportion to the section area of material 
removed (4). For this reason, in many practical applications of high- 
speed gas turbines and jet engines, solid disks without bores are adopted 
when the disks are made integral with the shaft. Of course, the magni- 
tude of o,,, which is directly proportional to the square of the speed, 
is also increased with the increasing speeds. 


? The boldface numbers in parentheses refer to the references appended to this paper. 
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Based on the Tresca criterion and its associated flow rule, and using 
the power function stress-creep-rate relation, the constant thickness solid 
disk at the steady-state condition has been treated by Wahl (5). In the 
present paper, on the basis of the same criterion and flow rule, we use 
the exponential function stress-creep-rate relation at steady-state con- 
dition. The variable-thickness solid disk will be analyzed and dis- 
cussed. Examples applying to the theoretical analysis for stress dis- 
tributions of the disk will also be given. 


2. CREEP RATE AND STRAIN EQUATIONS 


In general, for a given material, creep rate is primarily a function of 
stress, time and temperature, as given in the previous paper (1). Now, 
on this basis it follows that creep rate é is equal to the product of a 
function of stress, F(o), a function of time, f(t), and a function of tem- 
perature g(7), or 


e = F(a) f()g(7). (1) 


It is further assumed that the function of stress, F(o), may be repre- 
sented by 


(a) the power function creep law, o”, for relatively low stresses or 
low temperatures, 

(b) the exponential function creep law, e7/”, for relatively high 
stresses or high temperatures, 


where o is the equivalent or effective stress, m the ‘‘material creep 
number,” and a» the ‘‘material creep modulus.” 

Based on the basic assumptions given in the previous paper (1), and 
that the yield condition obeys the Tresca criterion and its associated 
flow rule, Eq. 1 may be written as follows: 


— é = kore!Tf(t), for > o2> (2a) 
= 0 


é = f(t) 
és = (1 —A)kore/7 f(t) for o:=02>0; and O<rA<1} (2) 
és = — kore-*!7 f(t) 


for the power function stress-creep-rate relation, and 
é3 = f(t), = 0 for > 62 > G3: (3a) 
é, = f(t) 


= (1 f(t) for o1 =o2>0,; and (3d) 
é; = — f(t) 


‘ 
oh 
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for the exponential. function stress-creep-rate relation, where \ is a 
multiplier or a parameter which will be further defined later, and the 
constants k and K may be determined from creep test data. Upon 
integration of Eqs. 2 and 3, the creep strains are, respectively, obtained. 


— € = kom f f(t, for (4a) 
0 


ko" f f(t)dt, etc. for 62> 0; (4d) 
0 


0 


etc. 
0 
for and O<A<1 = (5d) 


Fic. 1. 


in which the initial creep strains are assumed to be zero at ¢ = 0. 
Remember that for the steady-state condition, the stresses are inde- 
pendent of the time, while the temperature and the time may depend 
upon each other. 

As mentioned before, turbine disks of jet propulsion units are 


| 
and 
= 
O- + 
A 
\ 
Sid 


412 B. M. Ma [J. F. 


usually operated under the condition of relatively high stresses and high 
temperatures. For this reason in the present paper, we shall limit 
further discussion to the exponential function stress-creep-rate rela- 
tion only. 


3. CREEP ANALYSIS IN ROTATING DISK WITH VARIABLE THICKNESS AND 
UNIFORM TEMPERATURE—STEADY-STATE CONDITION 


On the basis of the basic assumptions given in the previous paper (1) 
the principal stresses in the disk, Fig. 1, may take the following biaxial 
stress system: 


¢,=¢ for 
Now, substituting this sytem into Eq. 3, we may have (1) 


Ke t —a/ f(t) 
TX 


for >o, and x.<x <1; 
u 


= AK otl/oop—a/T t 


Ou 


(1 — f(t) (10) 


= — (& + &) = f(t) 
for o.=0,, x1 <x <x and }<A<1 (11) 


in which \ is a multiplier and a function of x (or r). It is easy to see 
that on the basis of the incompressibility equation, that is é, + é- + €. = 0, 
h~t atx =x; ~Oandd\ = 1atx =x.. 

Consider now the element at radius 7 at the disk, Fig. 1, the equa- 
tion of equilibrium is given by 


d(hxa,)/dx — ho, + o,hx? = 0 (12) 


where / is the thickness of the disk, which varies with x (or 7), and 
o, = pw*r,?, the equivalent inertial stress in rotating thin ring of ra- 
dius 7». 


Outer Region of the Disk a, > 0, > 0, x. <x <1 


Differentiating Eq. 7 with respect to x, in which substituting from 
Eqs. 7 and 8, respectively, for w/r,x and d%/r,dx, it follows that 1/x 


z = C3. (6) 
4 
- 
Ou 
= = ( ( 8) 
r 
and 
‘ 
ag 
: 
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= —dao./odx. Integrating this between x, and x, we have the non- 
dimensional form for the tangential stress 


= — In (13) 


in which o;, is the value of ¢, at x = x... Now, consider the hyperbolic 
profile of the disk with the form h = cr,™x™ (where c and m are constant, 
but m is a negative value). By substituting this thickness function and 
Eq. 13 into Eq. 12 and integrating between x, and x, we have for the 
radial stress 


1 xltm + mx Ot. 


It may be readily verified that at x = x., Eqs. 13 and 14 reduce identi- 
cally to o,, = ¢,, which is to be determined. Here from the mathe- 


matical point of view the expression for the disk thickness 4 = cr,™x™ 
might be open to invite criticism. Because when x = 0, / will theo- 
retically become infinite. From the physical point of view, however, 
this explanation of the situation may be eventually used to interpret 
the fact that for a solid disk integrated with a shaft, the center of the 
disk is the shaft which may be considered as infinitely thick compared 
to the thickness of the disk proper. In the same manner, by integrating 
Eq. 12 between x = x. and x = 1 and rearranging, we obtain 


Oa 1 + mx,'+™ Ca 1 +m In + 1 m 3+m (15) 


which can be used to calculate o;,/¢. when the values of o,,/@4, 00/00, 
m, and x, are given. It may be noted that for the special case when 
x, = 1, Eq. 15 reduces to o:, = o,, = o,, at the outer radius of the disk, 
r,. Nowassuming m = — 3, and = 0.0275 (this is corresponding 
to the assumption that a» = 7,000 psi, and o, = 254,000 psi at the speed 
of the disk 17,800 rpm and the specific weight of the material y = 0.284 Ib 
per cu. in.), the values of o;,/0. versus x, are plotted for various values 
of a,,/¢,, Fig. 2. 

It is seen that in this figure the peaks of o,,/o, versus x, shift to the 
right or outward from the center of the disk as the values of a,,/a. 
increase. It may be mentioned here that form = — 3} andr = 10, the 


: 

=f. 
1 x xltm l+m : 
is c 

1+m Xe 1+m 

xitm — x dtm 

i = Xe < x < 1 ° ( 14) 

3+m 

ee : 


414 B. M. Ma [J. F. 1. 


thickness distribution of the disk profile is similar to that shown in 
Fig. 1. 

To determine the value of the critical point x, for which the condition, 
a, = a, > 0, is fulfilled and thereby Eqs. 13 and 14 will not hold any- 
more. Now, we divide Eq. 14 by Eq. 13, differentiate this ratio o,/o, 


Or o/ Og 0.5 


0 


0 02 04 O6 O8 10 
Uo 


Fic. 2. Tangential stresses at x, for various values of oro/oa, m = — 4, o0/¢a = 0.0275. 


with respect to x and set it equal to zero for a maximum at x = x. 
It yields 


(16) 


with given values of o:,/¢4, 00/00, Yo and m, x, can be readily found. 
For example, let o:,/¢. = 0.176, o0/o0. = 0.0275, and m = 

critical value of x. is approximately equal to 0.340. In this particular 
case, any value of x, less than 0.340 will not meet the required condi- 
tion, ¢, > o, > 0, over a portion of the outer region of the disk. 


05 
" 
= 
Oa MX Ga 
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Inner Region of the Disk o, = 6, > 0, S 1 
In this region of the disk, since 0, = o, for x; < x < x., Eq. 12 may 
be integrated between x and x... Upon substitution of h = cr.™x™ and 
integration, we have 
C1 Be m C1 
Ga Oa ( x ) Oa (2 + m)x™ ( 
In order to determine the value of \ if it satisfies the condition 
<1 forx; <x < at first, we differentiate Eq. 9 with respect 


0.6 


J 0.5 


Oro / Oa ® 0.2 


0.4 


0.8 1.0 


Fic. 3. Stress distribution in the disk m = — }, x. = 0.5. 


to x, (remember that \ is a function of x), next, we use Eqs. 9 and 10 
for substituting w/r.x and du/r.dx as before and finally we rearrange, 


getting 
dy/dx = (2/x + da./adx) = 1/x. 


Integrating this between x and x, and taking \ = 1 at x = x,, we have 


1 
+ ——— x; <x <x, (18) 


x 


in which o,, is given by Eq. 15 and o, is given by Eq. 17. By substi- 
tuting the value of o, from Eq. 17 into the integrand of Eq. 18 and 


| 

Fd 
X= 


6; /60 
0.4 


6,./6,= 0.2 


/ 
| 


02 O04 O86 1.0 
X= 


Fic. 4. Stress distributions in the disk m = — 3, x, = 0.6. 


0.7 


0 0.2 0.6 0.8 


Fic. 5. Stress distribution in flat disk m = 0, x. = 0.5, 
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introducing a new variable the integration for \ for given value of m 
may be found. Experience shows, however, that by a numerical 
evaluation of the integral in Eq. 18 it is relatively simple to determine 
whether the condition } <} <1 holds for x; <x <x. or not; a 
check reveals that this is usually true. 

The procedures for the creep analysis of variable-thickness and uni- 
form-temperature disk in the present paper may be summed up as 
follows. 


For the outer region of the disk, o, > o, > 0 forx, < x <1: 

(a) Use Eq. 15 to determine o;,/¢, for prescribed values of o;,,/oa, 

(b) Check on the value of o;,/0. from Eq. 16. 

(c) Substitute this value of o,,/0. into Eqs. 13 and 14, then the stress 
distributions of o,/o, and o,/o, are determined and the condition 
a, > a, holds. 


For the inner region of the disk, o, = 0, >0,x; < x < xX: 

(d) Use Eq. 17 to determine o,/o. = o,/o, after the value of o;,/o, has 
been substituted from Eq. 15. 

(e) Check on the value of \ from Eq. 18 if the condition 3} < \ < 1 for 
x; <x <x, is satisfied. 


From the above procedures, the stress distributions of the variable- 
thickness and uniform-temperature disk having m = — 3, r. = 10 in., 
oo/¢, = 0.0275 and temperature at 1200 F. for the various values of 
o,,/0, and x, are calculated and plotted in Figs. 3 and 4, respectively. 

For the purpose of comparison the stress distributions of a flat disk, 
m = ( (the other data remain the same), are plotted in Figs. 5 and 6, 
respectively. 


4. DISCUSSION AND CONCLUSIONS 


From the foregoing analysis, the following discussion and conclusions 
may be made. 


1. Based on the Tresca criterion, from Eq. 3a the effective creep 
rate é corresponding to the effective stress ¢ is given by 
= f(t). (19) 


At the steady-state condition and the given temperature, the stresses 
will be independent of time and temperature. Then, by differentiating 
Eq. 19 with respect to a, the slope of the creep curves for the exponential 
function stress-creep-rate relation is obtained. 


di/do = (20) 


It is seen that for a given material, by using Eq. 20, the value of creep 
material modulus a) may be readily determined from simple tension 


| 
; 
: 
| 
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creep-test data. In the previous example the value of oo is taken 
approximately for the disk material of Cr-Ni alloy steel at 1200 F. (6). 
2. From Figs. 3 and 4 for the variable-thickness disk of m = — }, 
it is seen that the peak stresses in the disk increase with the values of 
,,/¢, and x, and have a tendency to shift toward the outer radius 7,. 
3. Equations 14, 15, 16 and 17 for the variable-thickness disks can 


be immediately reduced to that for the flat disk when m = 0, Eqs. 13 


0.8 


0.7 


64/6, =6r/6, 
| 
6,./64=0.2 


0.6 


0.5 


0.4 


0.3 


0.2 


C2 0.4 0.6 . 1.0 
Xs 


Fic. 6. Stress distribution in flat disk m = 0, x, = 0.6. 


and 18 remain in the same form, but the numerical value of o:/¢. 
changes with that of o,,/0, according to the changes of Eqs. 15 and 17, 
respectively. The maximum stresses in the flat disk, Figs. 5 and 6, 
occurred at the center of the disk are as anticipated. They are quite, 
however, contrary to that of the variable-thickness disk, Figs. 3 and 4. 

4. There are two possible, extreme cases of the stresses distributions 
for variable-thickness disks (similarly for flat disks when m = 0 as 


4 
= » / 6a =0.4 
6+/ Ca 
5 
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discussed in the preceding paragraph) : If 


(a) o, > o, > Oat x; < x < 1, then only Eqs. 13 and 14 are needed 
to calculate the stresses. 

(b) o. =o, > 0 at x; < x < 1, then only Eq. 17 is required for 

determining the stresses. 


5. It is noticed in this creep analysis that based on the Tresca 
criterion and its associated flow rule, from the equilibrium equation, 
incompressibility equation and the stress-creep-rate relation the stress 
equations are derived in the simple closed form. 

6. Experience shows that the values of the maximum tangential 
stress based on the Tresca criterion and its associated flow rule, on one 
hand, are considerably higher than the test data while the values pre- 
dicted by the Mises criterion and its associated flow rule, on the other 
hand, are appreciably lower than the test data (2). It is felt that fur- 
ther spin tests of disks at elevated temperatures for various materials 
are needed before general conclusions can be drawn. 
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THE FRANKLIN INSTITUTE 


MINUTES OF THE STATED MEETING 
April 20, 1960 


The Stated Meeting of The Franklin Institute was held at 8:15 p.m. in the Lecture Hall. 
Mr. J. G. Richard Heckscher, Executive Vice President, called the meeting to order with 
approximately 160 members and guests present. 


Mr. Heckscher stated that the minutes of the February 17, 1960 Stated Meeting were 
published in the March issue of the JouRNAL. There being no corrections or additions these 
minutes were approved as published. The minutes of the March 16, 1960 Stated Meeting, 
he stated, would be published in the next issue of the JOURNAL and would be presented at a later 
meeting for approval. 


Before presenting the speaker of the evening, Mr. Heckscher noted that tonight we were 
presenting the Eleventh Charles Day Lecture, which is an endowed lectureship established in 
June 1941 providing for lectures on engineering and science. Mr. Day was a member of The 
Institute from 1900, becoming a life member in 1910. He was a member of the Board of 


Managers from 1908 until his death in 1931 and a Vice President of The Institute from 1920 
t» 1923, 


Mr. Heckscher then introduced Professor Stanley Keith Runcorn, Director of Physics 
Department, King’s College, Newcastle upon Tyne, England, whose subject was ‘Rock 
Magnetism, Polar Wandering and Continental Drift.” 


Professor Runcorn’s talk highlighted the most recent advances in the science of geophysics 
and in the course of his talk he used many slides of diagrams covering a broad phase of the sub- 
ject and excellent color slides showing important and historical rock formation. 


An interesting question period followed Professor Runcorn’s talk, and the meeting was 
adjourned by Mr. Heckscher at 9:45 p.M., after he had expressed on behalf of The Institute 
our very deep thanks to Professor Runcorn for his inspiring talk. 


WILLIAM F. JACKSON, JR. 
Secretary 
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BOOK REVIEWS 


ELECTROMECHANICAL ENERGY CONVERSION, 
by David C. White and Herbert H. Wood- 
son. 646 pages, illustrations, 6 X 9 in. 
New York, John Wiley & Sons, Inc., 1959. 
Price, $12.50. 


The generalization of the analysis of elec- 
tromechanical systems is a welcome feature 
of this book. In as much as we introduce the 
more intangible electromagnetic and electro- 
static field theory, in the analysis of reactions, 
couplings and energy transfer from electric 
circuits to mechanical systems, the statement 
of suitable parameters or coordinators in the 
dynamical approach is important. 

To be sure, the application of generalized 
coordinates to electromechanical systems and 
the use of Lagrange’s equations are well known 
in the literature, but the authors have credit- 
ably extended its applications to many in- 
teresting problems, which hithertofore have 
been approximated by piecemeal circuit 
analysis. 

In connection with the coupling reactions 
between electrical and mechanical systems, 
in place of the usual explicit introduction of 
mutual and self-inductance as coefficients 
describing the properties of the electromag- 
netic field, the authors have directly presented 
the parameters in terms of the flux linkages 
and the geometrical parameters and have ex- 
tended the analysis to nonlinear electromag- 
netic systems, where the energy /7(A, x)dA 
or the co-energy / (i, x)di are utilized. In 
its application for computing electromechani- 
cal forces, the authors have carefully presented 
the distinction between these different forms. 

In the authors’ “Fundamental Relations in 
Electromechanics,’’ the expression for D’Alem- 
bert’s principle is confusing, if not incorrect. 
The concept of D'Alembert’s principle is not 
postulated on the equilibrium of forces but 
rather on the equation of virtual work ex- 
tended to the inclusion of the inertia forces. 
Its extension to dynamical applications, as 
with Lagrange’s equation or in direct ten- 
sor form, is postulated on a choice of in- 
dependent coordinates consistent with the 
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constraints of the system. The authors’ 
equation 2 (psi — fii) = O evidently refers 
to a system of forces acting on the 
Ox; 


coordinate, where px = = 
x 


Now the generalized inertia force referred to 
Oxi 
the generalized coordinate x; is 2mz; _ and 
Xk 
mz; — }, since — is, in gen- 
OX’ OX’ 
eral, a function of the several independ- 
ent or generalized coordinates of the system. 
Thus D’Alembert’s equation should read 


Ox; Ox; 
where p, = Umi; —. Since Umz;— = pu, 
OX’ 


the authors’ equation Lz; = 0 seems mean- 

ingless. With nonholonomic systems the 

additional constraint reactions refer to con- 

dition relations 2 A,t, = 0, for the assumed 
(k) 


non-integrable coordinates defining the mo- 
tion of the system. 

Actually, the dynamical interpretation to 
closed circuits corresponds to generalized 
coordinates and the nature of the forces is 
quite different. Moreover, such coordinates 
are characterized as cyclic (or momentum) co- 


d oT 
Then in the equation aed 


ordinates. 
dt age 


= F,, the term ot should be deleted. The 


aV 
generalized force is Fy = — “9 +f. Then 
k 


we have » independent equations in the form 
pe — Fx = 0, for the n loops, where the cou- 
plings between the circuits may appear in any 
one of the terms, Pe or or explicitly in 
dt 
fii. This gives a physical interpretation and 
not an analogy. It seems, therefore, unnec- 
essary once the generalization is established, 
to introduce such terms as the inertia force 
of the kth node and the z, velocity of the kt 


: 
Ox dp, aT 
OX’ dt OX’ 
ie 
E 
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loop when circuit coordinates physically 
characterize the system. 

Extending T to include the complete elec- 
tromechanical system, then, since some of the 
coefficients of the electrical coordinates g, may 
be a function of the mechanical displacement 


oT 
x, then accounts for the electromechanical 
x 


forces. It should be pointed out that the 
motional emf’s, dependent on the mechanical 
velocities , which are time derivatives of the 
coefficients with respect to x, of the momen- 
tum terms of the electric circuits, account for 
the electromechanical energy transfer and are 
of a gyroscopic nature, with respect to the 
electromechanical forces. These couplings 
disappear in the over-all energy equation. In 
like manner, the electrostatic capacity coeffi- 
cients of g when a function of x then includes 


aV ‘ 
Th for the electrostatic-mechanical forces, 
along with the elastic mechanical reactions; 
aV A 
while — , correspond to the capacity emf’s 


for the electric circuits. 

The authors state the Lagrangian T— V=L 
as not sufficiently general to allow it to be 
extended directly to nonlinear systems. Ac- 
tually, Z applies to either linear or nonlinear 
systems, though not, without modification, 
to non-conservative systems. 

A major theme of this book includes a study 
of coupling forces and the energy transfer 
between the electric circuit system and the 
mechanical parts with special reference to non- 
linear relations of the parameters of the elec- 
trostatic or electromagnetic field systems. 
The electromechanical forces are expressed in 
terms of either the co-energy or energy form, 


where f, = — With lin- 


ear relations, these reduce to fj =i; 2 — 
r=1 


=f, Mee between pairs of circuits 7: and 7, with 
1 
mutual inductance M,, where Air = Miriy. 
The corresponding reactions on the electrical 
circuits are expressed in terms of the rate of 
change of the flux linkages. The field flux 
linkage is a function of both the current exci- 
tation and the mechanical displacement, 
where the latter is associated with the motional 
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counter emf and the corresponding electro- 
mechanical energy transfer. The mutual 
coupling between circuits renders this method 
of approach particularly applicable to a 
matrix tabulation, that is, a tensor analysis of 
electrical machines. 

The authors have extended this procedure 
to the analysis of various types of transducers 
and to generalized rotating machines, followed 
by specific applications to commutator ma- 
chines and to a-c. machines as induction 
machines. In conjunction with the tabula- 
tion for complex rotating machine circuit 
analysis, the authors include some discussions 
on system dynamics and control analysis and 
corresponding block diagrams, with the an- 
alytical techniques as solutions by Laplace 
transforms. 

In the authors’ Generalized Rotating Ma- 
chines a systematic analysis is developed 
through a well-defined grouping of windings. 
Two phase belts on stator and rotor with 
excitation are distributed sinusoidally at a 
electrical radians. The two phase belts are 
at quadrature. Arbitrary a* and b* stator 
magnetic excitation axes at quadrature are 
used as reference axes. The mechanical 
rotational angular coordinates # establishes 
the phase relation of the rotor quadrature coil 
beltsaand 6b. Obviously this can be extended 
to any number of phases ultimately reduced 
to quadrature phases. To arrive at the self- 
and mutual inductance coefficients between 
coils, it is necessary to establish the air gap 
flux distribution. This is calculated from the 
magnetic energy Wa = 4/HBdv in terms 
of the circuit parameters 422L;;i;7; where 
special consideration is given to pole saliency, 

From this rotating machine configuration. 
the four circuit voltage equations, that is for 
the quadrature a and b phases of the rotor and 
stator, respectively, and the mutual torque 
in terms of the rate of change of the mutual 
flux linkages of the several circuits with 
respect to the rotor mechanical displacement 
@. An interesting feature for the derivation 
of the latter is the use of the derivative of the 
over-all work equation, where the magnetic 
energy Wy is first stated in terms of the 
energy form and is reduced to its co-energy 
form Wy’ = SXA;(i, ¢)di;. The method is 
also explicitly stated in its Lagrangian form, 
where the circuit and_ electromechanical 
coupling analysis are considered as a system. 
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Certain physical relations of the induction 
motor are compared with a fluid coupling. 
Actually, the relations are quite close. The 
working mediums are respectively the air gap 
flux or the fluid, which transfers the mutual 
torque and rotor. The 
slippage energy Ts-w* goes into heat through 
the /?R losses in the rotor or in the heating 
of the fluid with the clutch. The torque reac- 
tion times the synchronous rotating field is 
equivalent to the counter emf times the stator 
currents, or to the authors’ stator power 
(Eq. 3-239); and, with the clutch, this is its 
mechanical input from the drive shaft. These 
relations also concur with the authors’ Eq. 
3-237, where part of the slippage energy is 
transferred externally from the rotor circuit 
in the term V’J" cos 6’. The potentiality of 
the induction motor as a “‘controlled”’ slippage 
clutch is great. 

Undoubtedly the approach through tensor 
formulation of the interactions and energy 
transfer between complicated circuit systems 
is particularly adapted for computer methods, 
but such methods are limited in the over-all 
interpretation of the physical characteristics 
of electrical machinery. The concept for 
over-all approximations of field excitation, 
flux distribution, torque and power relation, 
in terms of distributed current excitation with 
stationary and revolving air gap waves; and 
the space phase relations of corresponding 
flux distributions, with ampere turns along 
the gap, in conjunction with various types of 
windings over-all macroscopic 
parameters; permit equally important gen- 
eralizations for various classes of machines 
with a physical picture of the interactions. 
This approach is particularly useful in a-c. 
machinery as with the induction motor. It is 
important to distinguish the air gap torque 
distribution and relations of space and time 
phase of the flux and counter emf's. The 
relations of various methods of excitation, 
armature reactions, with the above, are neces- 
sary for analytical design. The coordinating 
of these methods with the authors’ tensor 
analysis would be an interesting undertaking. 

It is to be noted the circuit equations and 
the electromechanical coupling forces can be 
established directly for the majority of 
examples listed. It is equally important to 
parallel a clear understanding of the nature of 
the reactions and energy transfer between 


between stator 


treated as 
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circuits though a step-by-step build-up of 
the interacting electromagnetic forces and to 
distinguish directly the motional and trans- 
form emf’s relating to the mechanical energy 
transfer and that to the field. Further, 
consideration should be given to the field 
excitation and flux build-up for the magnetic 
circuit and the gap distribution in relation to 
the electromechanical forces and hysteresis 
torque. This latter phase is one of the most 
difficult problems in the design of electro- 
mechanical systems. 

From the reviewer's point of view, it would 
seem a more simplified approach of the simple 
basic characteristics, as to the nature of 
torque, current and flux, motional and trans- 
former counter emf’s, field excitation relations, 
rotating fields, etc., should precede the more 
ambitious program in the use of Lagrange’s 
equations and such matters as quasi coordi- 
nates. The authors should note a large 
bibliography on the subject which deserves a 
careful scientific review. 

The above is not suggested as a criticism of 
the authors’ work, which enjoys an important 
contribution to the art, but rather as a criti- 
cism for a curriculum in teaching this subject. 

Moreover, this book deserves special com- 
mendation in the development of a systematic 
and most interesting generalized analysis of 
electromechanical energy conversion systems. 


RupEN EKSERGIAN 
The Franklin Institute Laboratories 


MECHANISMS AND Motion, by K. H. Hunt. 
114 pages, diagrams, 54 X 8} in. New 
York, John Wiley & Sons, Inc., 1959. 
Price, $4.25. 


Since the subject of kinematics is an indis- 
pensible preliminary of such fields as dynam- 
ics of machinery and plastic deformation of 
frameworks, it is unfortunate that it should 
be so poorly treated by writers of engineering 
textbooks. Most such books contain such a 
mixture of vectorial dynamics, drafting room 
exercise, practical design limits and illustra- 
tions of mechanisms of varying complexity 
that the portion devoted to kinematics is 
often barely sufficient to justify the title. 

While the writer of the present volume has 
made a laudable attempt (especially in the 
last few chapters) to present the abstract 
aspects of kinematics, from the point of view 
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of the machine designer, his treatment of 
fundamentals can hardly be recommended. 
His discussion of so basic and important a 
concept as degrees of freedom, to which 
several pages are devoted, is quite hazy com- 
pared for example, to the brief and intuitively 
clear exposition in Sommerfeld’s ‘“‘Mechan- 
ics."’ Vector methods are discussed in one 
brief chapter which contains the erroneous 
statement that in plane motion, the relative 
centripetal and Cariolis components of ac- 
celeration are parallel, which is not true if the 
point in question has a radial component of 
velocity in the rotating coordinate system. 
The elaborate vector methods are hardly 
justified in the solution of extremely single 
examples in which the author fails to make 
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clear an aspect of the Cariolis acceleration 
which e.g. Webster’s treatment brings out; 
i.e., that it appears as the cross product in the 
expansion of the square of a sum. The au- 
thor’s solutions are thus quite unnecessarily 
complex especially in the case when the axes 
of rotation of the fixed and moving systems 
of coordinate are coincident. 

The final chapters on space mechanics and 
advanced kinematics of plane mechanisms 
are more worthwhile, although the use of 
two-dimensional symbols for hinged guides 
in three-dimensional pictures does not con- 
tribute to clarity. A number of interesting 
mechanisms are shown in the problems. 


CAPPEL 
The Franklin Institute Laboratories 
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MATERIALS AND TECHNIQUES FOR ELECTRON TuBEs, by Walter H. Kohl. 
New York, Reinhold Publishing Corp., 1960. 


638 pages, diagrams, 6 X 9 in. 
$16.50. 


Fa._out, edited by John M. Fowler, with a foreword by Adlai E. Stevenson. 
New York, Basic Books, Inc., 1960. 


diagrams, 6 X 9 in. 


INDUSTRIAL INSTRUMENT SERVICING HANDBOOK, by Grady C. Carroll. 
New York, McGraw-Hill Book Co., Inc., 1960. 


tions, 6 X 9 in. 


ELECTRONIC COMPUTERS. 
263 pages, illustrations, 5} & 8} in. 
$15.00. 


PRINCIPLES AND APPLICATIONS, by T. E. IVALL. 
New York, Philosophical Library, 1960. 


Revised edition, 


Price, 
235 pages, 
Price, $5.50. 


848 pages, illustra- 
Price, $16.00. 


Second edition, 
Price, 


AERONAUTICAL DicTIONARY, compiled by the National Aeronautics and Space Administration. 


197 pages, diagrams, 5} X 9 in. 
$1.75. 


Washington, Government Printing Office, 1960. 


Price, 


GERMAN-ENGLISH MATHEMATICS DICTIONARY, compiled and edited by Charles Hyman. 


131 pages, 8} X 11 in. 
No price given. 


MATHEMATICS REFRESHER, by Kurt Wolter. 
Price, $3.75. 


Philosophical Library, 1959. 


ELEMENTARY MOoperN Puysics, by Richard T. Weidner and Robert L. Sells. 
Boston, Allyn and Bacon, Inc., 1960. 


illustrations, 6 X 9 in. 


Data Book For Civi_ ENGINEERS. 
v.p., illustrations, 9 X 11} in. 


96 pages, diagrams, 4} X 7 in. 


DesiGn, by Elwyn E. Seelye. 
New York, John Wiley & Sons, Inc., 1960. Price, $24.00. 


New York, Interlanguage Dictionaries Publishing Corp., 1960. 


New York, 


513 pages, 
Price, $8.50. 


Third edition, Volume 1, 
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VACUUM PROCESSING IN METALWORKING, by 
J. Wesley Cable. 202 pages, illustrations, 
5 X 74 in. New York, Reinhold Publish- 
ing Corp., 1960. Price, $5.50. 


In this small text the relatively new field 
of vacuum technology receives complete and 
practical treatment. It also encompasses a 
wide range of operations in the metalworking 
field. 
nology such as the means of producing low 
pressures and their measurement is reviewed. 
The author's present vacuum 
technology in a descriptive manner, rather 
than in a highly technical form. As a result, 
it is beneficial to any industry involved in the 
refining and processing of metals. 


The general phases of vacuum tech- 


aim is to 


ROCKET PROPELLANT HANDBOOK, by Boris 
Kit and Douglas S. Evered. 354 pages, 
diagrams, 6} X 9} in. New York, The 
Macmillan Company, 1960. Price, $12.50. 


Nearly a hundred chemicals having possible 
use as rocket propellants are described in this 
concise handbook. The chemicals include 
fuels, oxidizers, and monopropellants. Tables 
are included giving such data as heat value, 
density at various temperatures, specific im- 
pulse and performance with alternative oxi- 
dizers or fuels. The text will be of interest 
not only to the rocket propulsion industry 
and the military services, but also to the 
chemical industry, manufacturers of handling 
and storage equipment, and engineering com- 
panies involved in the design and construction 
of facilities. 


AIRCRAFT AND MISSILE DESIGN AND MAIN- 
TENANCE HANDBOOK, by Charles A. Over- 
by. 369 pages, diagrams, 6 X 9 in. New 
York, The Macmillan Company, 1960. 
Price, $9.75. 


A practical handbook for mechanics, in- 
spectors, draftsmen, supervisors, designers 


and modification engineers. Containing over 
160 tables and 40 figures to illustrate the 
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standard methods of equipment installation 
and maintenance, the book gives standards 
and hundreds of tips on the proper handling 
of materials and tools. Topics covered are 
electrical and plumbing systems, materials of 
construction, aircraft and missile hardware, 
color codes and conversion systems, and 
processes such as metal spraying, anodizing, 
rust-proofing, plating, welding and brazing. 


ELEMENTS OF ION EXCHANGE, by Robert 
Kunin. 164 pages, diagrams, 5 X 7} in. 
New York, Reinhold Publishing Corp., 
1960. Price, $5.75. 


This little book provides a basic under- 
standing of ion exchange techniques, written 
specifically for those with little or no theo- 
retical background. As there is scarcely an 
industry today in which ion exchange is not 
involved, this concise handbook, giving the 
description of its development and importance, 
will be of valuable assistance to executives 
and operating personnel. In addition, a broad 
view of the subject is provided for analytical, 
physical and industrial chemists, biochemists, 
agronomists and nuclear researchers. 


THE AERODYNAMICS OF POWERED FLIGHT, by 
Robert L. Carroll. 275 pages, diagrams, 
5} X 9in. New York, John Wiley & Sons, 
Inc., 1960. Price, $8.50. 


The author gives a new approach to the 
subject, both simple and direct. By showing 
that lift is directly dependent on power and 
by treating rocket power as the logical lift 
source, the book provides an exposition of the 
science and mechanisms of powered flight 
that is unique in its clarity and simplicity. 
A foundation in fluid flow and airfoil theory 
is detailed, while also introducing the funda- 
mentals of propulsion, stability, performance, 
design and control. In order to emphasize 
the practical aspects of flight, 216 problems 
are included. 
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PETROLEUM ENGINEERING: DRILLING AND 
WELL CompLetions, by Carl Gatlin. 341 
pages, diagrams, 84 X 11 in. New York, 
Prentice-Hall, Inc., 1960. Price, $13.00 
(trade edition); $7.95 (text edition for 
classroom adoption.) 


The technical as well as non-technical reader 


will find this book useful as a reference and > 


practical for self-study of the problems of 
drilling and completing oil wells. The author 
examines the basic elements of drilling, test- 
ing, formation, evaluation and completions, 
also surveying various geological and geo- 
physical exploration methods. An extensive 
bibliography, listed by chapters is an excellent 
guide to special topics. Illustrations com- 
bined with an uncluttered text provide a 
clear, readable coverage of well location and 
development. 


Acry.ic Resins, by Milton B. Horn. 184 
pages, illustrations, 5 X 7} in. New York, 
Reinhold Publishing Corp., 1960. Price, 
$4.50. 


This little handbook in the Plastics Appli- 
cations Series, describes the four types of 
acrylics in respect to their manufacture, fabri- 
cation and applications. These types include 
cast products, molding compounds, emulsion 
and solution compounds. The author explains 
raw materials and manufacture of these 
polymers, then he covers the current and 
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future developments in terms of potential 
use in a great variety of products. Thus, a 
complete guide is given to chemists, sales 
personnel and executives on every aspect of 
acrylics in addition to a broad understanding 
of uses and future trends. 


Puysics: FoR STUDENTS OF SCIENCE AND 
ENGINEERING, by Robert Resnick and 
David Halliday. 554 pages, diagrams, 
6 X 9} in. New York, John Wiley & Sons, 
Inc., 1960. Price, $6.00. 


The authors have aims to reveal the unity 
of physics in many ways throughout this text. 
For example, the conservation laws of energy, 
linear momentum, angular momentum, and 
charge are used repeatedly. Wave concepts 
and properties of vibrating systems, such as 
resonance, are used in mechanics, sound, 
electromagnetism, optics, atomic physics and 
nuclear physics. The field concept is applied 
to gravitation, fluid flow, electromagnetism 
and nuclear physics. The mathematical level 
of the book assumes a concurrent course in 
calculus. An unusually large number of 
worked-out examples, with the ‘plug-in’ 
variety used only to emphasize a numerical 
magnitude. Algebraic, rather than numerical 
solutions are stressed. The mks system of 
units has been adopted throughout the text, 
although the British engineering system is 
also used in mechanics. 
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Tiny Microphones Record Space 
‘Hits’. The 142-lb. Explorer VI satel- 
lite is carrying two tiny microphones 
as part of a micrometeorite detection 
experiment sponsored by the Air Force 
Cambridge Research Center (AFCRC), 
Bedford, Mass. AFCRC is an Air 
Research and Development Command 
(ARDC) facility. 

Explorer VI, launched August 7, 
1959, circles the Earth at an estimated 
22,000 mph. Its orbit was_ pro- 
grammed for 23,000 statute miles at 
apogee (farthest point from Earth) 
and 160 statute miles at perigee (closest 
point to Earth). Scientists estimate 
the satellite’s life span at two years. 

Purpose of the instruments is to 
measure the number of times the 


“paddle wheel”’ satellite, launched un- 


der the direction of the National Aero- 
nautics and Space Administration 
(NASA), is struck by these tiny par- 
ticles in space. 

The microphones, located opposite 
each other, are firmly mounted in the 
center of 6-in. by 12-in. strips of the 
satellite’s ;'g-in. thick aluminum skin. 
They detect either a strong or weak 
hit by particles one micron (a thou- 
sandth of a millimeter) in diameter or 
larger. More than 25,400 microns 
placed side by side would measure but 
a single inch. 

One of these dust particles striking 
the satellite’s surface generates an 
acoustical pulse which, after traveling 
through the exposed plate, is trans- 
formed by the microphone into an 
electrical pulse. This electrical pulse 
is in turn transmitted to data accumu- 
lating stations on Earth. Size of the 
particles is determined between cer- 
tain known limits from the magnitude 
of the pulses generated. 
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Thermocouple Measures Tempera- 
tures of Molten Metal.—A new ex- 
pendable immersion thermocouple for 
measuring molten metal temperatures 
has been developed by the Instru- 
ments and Systems Division of Engel- 
hard Industries, Inc., Newark, N. J. 
Primarily designed for measuring steel 
temperatures in both open hearth and 
electric furnaces, it may be used with 
equal efficiency with other molten 
metals. The new unit combines im- 
proved economy, faster response, high 
accuracy and simplified assembly. 

The expendable thermocouple de- 
sign is simple consisting of a ceramic 
base which serves as a mounting for 
the noble metal thermocouple wires. 
A firmly anchored transparent fused 
silica tube protects the thermocouple 
wires and is, in turn, protected from 
slag and floating solids by a special 
metal cap. The cap absorbs the phys- 
ical shock of immersion and melts 
away immediately after it has plunged 
beneath the slag. 

Special design of the thermocouple 
unit and its receptacle allows the unit 
to be assembled without requiring 
alignment of contact pins. Design 
also assures that molten metal cannot 
penetrate the assembly and damage 
the plug receptacle. 

Two different thermocouple com- 
binations are available in the Engel- 
hard unit, either platinum vs. platinum 
10 per cent rhodium or platinum vs. 
platinum 13 per cent rhodium. 

Standard pipe, handle and lead wire 
assemblies may be used with the unit 
or these components may also be 
obtained from Engelhard. 
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Glass Fibers Used in Concrete.— 
A method of producing concrete struc- 
tural elements both reinforced and 
prestressed by glass fibers has been 
developed by the U. S. Army Engineer 
Research and Development Labora- 
tories, Fort Belvoir, Va. The method, 
which has been patented, was devel- 
oped by Mr. Solomon Goldfein, a 
chemical engineer, who is chief of the 
Plastics Section in the Materials 
Branch of the Laboratories. 

It is expected to be used in special 
applications, where the light weight, 
corrosion resistance, absence of mag- 
netic field, and electrical resistance 
of the reinforcement are important 
factors to be considered. It has been 
considered for such structures as pil- 
ing for docks exposed to salt water, 
concrete ships and vehicles requir- 
ing non-magnetic reinforcement, struc- 
tures where the reinforcement must 
not be allowed to conduct electricity, 
and as reinforcement for large concrete 
pipes conducting water. 

One of the major difficulties ex- 
perienced in producing such structural 
elements has been the problem of 
securing adequate bond of the glass 
fibers to concrete. The method de- 
veloped by the Laboratories solves 
this problem through the use of a 
binder, such as an epoxy resin, spe- 
cially formulated, which bonds the 
glass fibers to themselves and to the 
concrete, thus protecting them against 
the high alkalinity of the cement. 

The reinforced concrete may be 
made by several different methods. 
In one method, a prefabricated plastic 
glass fiber rod having a silicious or 
cementitious covering on its surface, 
is used. The rod is embedded in con- 
crete exactly as its steel reinforcement 
counterpart. The covering acts to 
increase the strength of the bond be- 
tween the glass fibers and concrete. 
In the other method, glass fiber im- 


CuRRENT Topics 


429 


pregnated plastic embedded in the 
concrete hardens during the setting 
of the concrete, thus effecting a bond. 

Prestressed and poststressed con- 
crete may be made by using principles 
established previously during the de- 
velopment of prestressed plastics by 
Goldfein. The rods may be prefabri- 
cated with special gripping devices 
bonded to them for tensioning. 

The material probably will never 
compete in ordinary building con- 
struction, because of its low modulus 
and high cost compared to steel. 
The former disadvantage, however, 
is expected to be considerably reduced 
soon by the use of a special glass 
fiber, which has a modulus almost 
equal to that of steel. 

In countries having no iron re- 
sources but plenty of glass-making 
raw material, the price differential 
will lose some of its significance. 


Leveling Agent Improves Dyeing.— 
Deeper, level and solidly dyed shades 
on wool and/or nylon stocks, yarns or 
fabrics can now be obtained econom- 
ically with neutral dyes by using a 
new leveling agent developed by 
American Cyanamid Company. 

Known as Calco Leveling Salt LV, 
the agent is applied with the dye 
under processing conditions, pH 4-6, 
that should be used for neutral dyes. 
The soft hand of wool is retained 
and the striped or Barré effect en- 
countered in dyeing nylon carpeting 
is eliminated. This salt is also ex- 
pected to solve similar problems in 
dyeing such other forms of nylon as 
hosiery and sweaters. 

Unlike other leveling agents, Calco 
Leveling Salt does not act as a re- 
tarding agent, which makes incom- 
plete exhaustion of the dye bath the 
penalty of satisfactory appearance. 
Cyanamid’s new agent allows the 
colorant to dye rather than stain the 
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cloth. As a result, shades may be 
matched by adding color in process, a 
procedure not possible with other 
agents of this type. 

Calco Leveling Salt LV serves the 
dual purpose of wetting out the fiber 
and controlling the dyeing. Its use 
does not have any adverse effect on 
the fastness properties of the dyes used. 


Improved Rubber Gaskets.— Ameri- 
can Felt Company, Glenville, Conn., 
has announced the introduction of 
“Vistex III,” a new and improved 
fiber-reinforced nitrile rubber gasket 
material. The new material, made of 
synthetic fiber-reinforced ‘Teflon,’ 
is recommended for its excellent 
stability to corrosive gases and liquids 
under both high and low temperature 
conditions. 

Among features claimed for the 
new ‘‘Vistex’’ are its 100 per cent 
compression recovery, its exceptional 
tensile strength, its resistance to oil 
and aromatic fuels, good flexibility, 
and its non-adhesive and non-corro- 
sive properties. 

“Vistex III’’ is available in sheet 
form, up to 30 in. X 60 in., and in 
standard 7s, # and thick- 
nesses. It has 84 lb. per cubic ft. 
density anda minimum tensile strength 
of 1000 psi. At 300° F., it retains 620 
psi. tensile strength and at 430°, loses 
only 105 psi. It is especially well 
suited for gaskets where light weight 
is a factor. One square yard of #5 
“Vistex”’ weighs only 30 oz. 

The new material can be preci- 
sion cut or punched into gaskets, 
strips or washers. It has excellent 
chemical stability to acids, alkalis, 
aliphatic and aromatic compounds. 
Among operating characteristics found 
to be inherent in the new material 
are excellent die and hand-cutting 
properties, freedom from fraying and 
raveling, and minimum-pressure seal- 
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ing on either polished or irregular 
surfaces. 

Under most conditions, it can be 
re-used. It has very high tear value 
and because of its high strength, it 
is an ideal industrial gasket material. 

Because of its high degree of chem- 
ical stability, and the versatility of 
uses, it is believed AFCO’s new “‘Vis- 
tex III’ will find wide application in 
a variety of industrial uses. 


The Ferreed.—A\ team of engineers 
at Bell Laboratories has invented a 
very fast electromechanical switch 
called the ‘‘ferreed.”’ It utilizes a new 
technique for actuation which is com- 
patible in speed with electronic circuits. 

The ferreed may find use as an inter- 
connecting element in telephone switch- 
ing networks, where it could be con- 
trolled a thousand times more rapidly 
than switches presently employed. 

Inventors of the new switch, A. 
Feiner, C. A. Lovell, T. N. Lowry, 
and P. G. Ridinger, are members of 
the technical staff of Bell Telephone 
Laboratories. They anticipate appli- 
cations for the ferreed in various 
types of systems where advantage 
may be taken of metallic contacts 
controllable by current pulses of a 
few microseconds duration. 

Magnetic materials that can be 
quickly switched between two alter- 
nate states have been widely used as 
‘‘memories” for storing information 
in digital computers. The useful out- 
put from these elements has generally 
been limited, however, to electrical 
signals of a transistory character. 

In conventional electromagnetic re- 
lays, on the other hand, continuous 
electrical currents can be used to open 
and close metallic contacts for ex- 
tended periods of time, but at speeds 
limited by the mechanical motion of 
the moving parts. 

The ferreed combines the rapid 
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switching of bistable magnetic ma- 
terial with metallic contacts for out- 
put indications that persist as long as 
desired without further application of 
power. In several models of the de- 
vice, a cobalt ferrite has been used 
as the magnetic material and a glass- 
sealed magnetic reed switch for the 
output contacts—hence the name, 
““‘ferreed.”’ 

In operation, the magnetic material 
is switched by a magnetomotive force 
applied, typically, as a five-ampere 
current pulse in a thirty-turn winding. 
Control pulses as short as five micro- 
seconds will switch the magnetic 
material, resulting in the passage of 
magnetic flux through the movable 
members of the reed switch. Actual 
closure of the contacts is delayed by 
inertia of the reeds for several hun- 
dred microseconds. 

Release of the contacts is brought 
about by cancellation of the magnetic 
flux through the reeds, as the result 
of another five-microsecond switching 
operation. Opening of the contacts 
requires less time than the closing 
operation. 

Several new magnetic materials 
have been synthesized for use in de- 
velopmental models of the ferreed. 
Among these are ferrites exhibiting 
characteristics midway between per- 
manent magnet* and computer mem- 
ory materials, and ferrite suspensions 
in plastic “tailor-made” for specific 
magnetic, electrical, and mechanical 
properties. The reed switch is similar 
to a type used in standard electro- 
magnetic relays for the past decade. 

In a telephone switching network, 
subscribers are interconnected by the 
closure, inspecified patterns, of switches 
called ‘‘crosspoints.”” A typical cen- 
tral telephone office contains many 
thousands of these switches; of these, 
a few dozen are operated for each call. 
Those crosspoints required for a par- 
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ticular connection must be selected 
by the coincidence of pairs of control 
pulses. The ferreed’s magnetic struc- 
ture responds to a pair of pulses by 
closing the output contacts. They 
are opened by a subsequent single 
pulse. 

Besides coincident pulse operation 
and rapid actuation, the magnetic 
structure permits the desired switch- 
ing by current pulses of widely vari- 
able character. 

The result, when ferreeds are used 
in switching networks, is_ reliable 
operation from simple control circuits. 


Fifteen-Story Missile Tower.—A 
15-story high portable missile service 
structure, believed to be the tallest 
and heaviest pneumatic tire-mounted 
structure in this country, has been 
developed for use at missile test sites, 
according to a recent Department of 
the Army announcement. 

The portable service tower was 
developed by the U. S. Army Engi- 
neer Research and Development Lab- 
oratories, Fort Belvoir, Va., for use in 
servicing missiles of the Redstone and 
Jupiter classes, or any missile up to 
136 ft. tall. The tower was designed 
and built by the Noble Company, 
Oakland, Calif. 

Previous missile towers have been 
either fixed or rail-mounted. The 
new 151-ft., 350,000-lb. structure, 
however, is mounted on two base 
trailers of six wheels each. The out- 
board rearmost wheels are driven by 
an electric motor, and front wheels are 
steered by hydraulically operated cyl- 
inders controlled by an electric sensing 
device running on a steering track. 

After being positioned, the two 
crane hooks operating from the “ham- 
merhead”’ top pick up the missile and 
set it on its launcher. Jack pads are 
lowered and the structure is off tires 
until the time for roll back. 
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The tower has six adjustable plat- 
forms for use by engineers and 
scientists in preparing a missile for 
flight. It also has two elevators, 
a complete intercom system, and a 
“panic button.” The button is de- 
signed for use at the most critical 
moment of missile fueling. If a fire 
should start, a push of the button, 
either on the tower or in the control 
bunker, slides the platforms back 
from around the missile, lifts up the 
jacks off the ground and wheels the 
tower itself away from the fire. 

An engine-generator, motor-genera- 
tor, hoists, motors, gear boxes, and 
electric control panels are mounted 
on the base trailers. The twin masts 


can be taken apart and the flat panels 
and platforms stacked on transport 
by aircraft, 


trailers for movement 
or for highway tow. 


Solid State Data Processing Sys- 
tem.—The solid state IBM 7080 data 
processing system has been announced 
by International Business Machines 
Corporation. It is the most powerful 
computer yet designed specifically for 
business with processing speeds up to 
ten times faster than the IBM 705 
systems. 

The 7080 is capable of producing 
management reports and processing 
payrolls, inventory control, billing 
and similar accounting jobs at a sub- 
stantially lower job cost than the 
previous systems. 

New design advances that give the 
7080 its high capacity include: (1) 
fastest memory available, called Com- 
munication Storage, which can 
transfer information between tape 
units and main data storage so fast 
(1.09 microsecond character rate) that 
up to five tape units can be reading 
in or out simultaneously; (2) a main 
memory from which any of up to 
160,000 characters of information can 
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be called into the central processing 
unit in 2.18 microseconds giving the 
7080 the fastest main storage avail- 
able; and (3) priority processing, a 
system which achieves maximum 
efficiency of simultaneous reading, 
writing and processing operations, 
permits input/output devices to auto- 
matically control the flow of programs. 

A typical 7080 system will sell for 
$2,528,000 or rent for $55,500 a 
month. This includes the central proc- 
essing unit, magnetic core storage, 
console, two tape control units and 
twelve magnetic tape units. A basic 
system will sell for $2,233,000 or rent 
for $49,100 a month. 

The new concepts in the 1.09 micro- 
second communication storage sim- 
plify the attachment of input/output 
devices and provide the basis for the 
priority processing system. 

As a result, the 7080 can take in or 
feed out information to or from 
sources such as its magnetic tape 
units, process and transmit data in- 
ternally, all at the same time. Com- 
bined with priority processing, this 
greatly reduces job times. | 

In one second, the 7080 can: read 
or write 312,500 characters of infor- 
mation; make 303,000 logical deci- 
sions; add or subtract 78,000 six-digit 
numbers; multiply 7,100 six-digit 
numbers. 

The 7080 has a unique, newly- 
designed console which includes a 
digital display of registers. This 
information is projected in back- 
lighted Arabic numbers for easy 
reading by the operator. Previously, 
this information appeared on systems’ 
consoles in the form of binary-coded 
lights. 

The 7080's increased main memory 
speed of 2.18 microseconds compares 
to seventeen microseconds on the 
IBM 705 I and II systems, and nine 
microseconds on the 705 III. A 
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magnetic core memory of either 
40,000, 80,000 or 160,000 characters 
is available on the new system. The 
705 I has a storage of 20,000 charac- 
ters, the 705 II has 40,000 and the 
705 III has an 80,000 character 
storage. 

Ultra high-speed access to memory 
in the 7080 significantly contributes 
to stepped-up computing. The 7080 
is designed to accept 705 I and II 
programs and internally process them 
ten times faster, and 705 III programs 
six times faster. Even greater speeds 
can be realized when 705 programs are 
modified to take advantage of the 
increased power and flexibility of the 
7080’s new program instructions. 

The new system has three modes of 
operation—705 I and II, 705 III and 
7080. The operator can establish 
either 705 mode simply by pushing 
a compatibility switch on the console. 

Major 705 programs which will 
operate on the 7080 system include 
Autocoder III, Input/Output Pack- 
age, Decision Making Language, 
FORTRAN, Report Generator, File 
Maintenance System, 705 Processor, 
Utility Programs and Sort and Merge 
Routines. 

The solid state 7080 needs only fifty 
per cent of the air conditioning and 
power of its predecessors, and thirty 
per cent less space. These reductions 
will permit expansion of a solid-state 
system within present facilities as 
capacity requirements grow. 

Even greater savings can be at- 
tained if peripheral tube equipment 
is replaced with the solid state IBM 
1401 system. Complete tape and 
tape unit compatibility is achieved 
between on-line and off-line equip- 
ment because the 729 II and 729 IV 
tape units are used in both systems. 

With 1401 off-line printing, punch- 
ing and card-to-tape conversion equip- 
ment, the 7080 becomes a powerful 
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system for applications such as de- 
mand deposit accounting and other 
banking operations; railroad freight 
car accounting; public utility custo- 
mer accounting and merchandise con- 
trol in the retail industry. 


Fluoroscopic Screen Detects Pres- 
ence of Stray Radiation.—A new-type 
fluoroscopic screen called RAD-A- 
LERT which enables radiologists to 
detect unabsorbed stray radiation, has 
been announced by United States 
Radium Corporation. A patent has 
The new screen— 
which provides this extra measure of 
safety and service to the radiologists— 
features a special border, a portion of 
which fluoresces red, in contrast to the 
yellow-green of the screen proper. 
This red signal indicates to the radiol- 
ogist that radiation is escaping be- 
yond the protective lead glass of the 
fluoroscopic assembly. 

The radiologist can immediately 
take steps to remedy this condition 
through reduction of shutter openings, 
altering _patient-screen _ distance, 
proper tube centering or other correc- 
tive measures. 

In actual use, the border of the 
RAD-A-LERT screen presents a level 
of brightness readily visible, yet not 
sufficiently bright to disturb dark 
adaptation or interfere with con- 
tinuous fluoroscopy. 


High-Speed Silicon pnpn Diode 
Switch.—A silicon pnpn diode capable 
of switching at extremely high speeds 
was described in a paper presented at 
the Winter General Meeting of the 
American Institute of Electrical Engi- 
neers. The switch, developed at Bell 
Telephone Laboratories, was a result 
of the joint efforts of A. N. Baker, 
J. M. Goldey, and I. M. Mackintosh. 

This device is useful in switching 
moderate amounts of power with 
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extreme rapidity. Typical applica- 
tions are in line pulsers, and high- 
speed gates for use in certain types of 
telephone equipment now under de- 
velopment. Turn-on time is approxi- 
mately ten millimicroseconds and 
turn-off time is about the same, in- 
cluding the storage time of the charge 
carriers—which is around 4 millimi- 
croseconds. More refined measure- 
ments show that current is down four 
orders of magnitude in less than 100 
millimicroseconds. Several watts of 
power can be switched by this device. 

The diode has a triple diffused 
structure, starting with a wafer of 
n-type silicon. The wafer is then 
etched to form a mesa-type construc- 
tion to reduce the area, resulting in a 
unit having a capacity of about 5 
micromicrofarads. 

High-speed operation is achieved 
by the use of thin base layers and by 
reducing the lifetime of the charge 
carriers. The latter can be done very 
precisely either by irradiating the 
device with high-energy electrons or 
by diffusing in minute amounts of 
gold. 

In the ‘‘on’’ condition, this diode 
has a resistance of about 2 ohms. Im- 
pedance when ‘‘off’’ is represented by 
a capacitance of 5 micromicrofarads 
with an extremely low leakage current 
of around a millimicroampere at room 
temperature. 


Undersea Arctic Cable.—An im- 
portant link in the rearward communi- 
cations for America’s Ballistic Missile 
Early Warning System (BMEWS) 
became a reality recently when the 
first undersea telephone cable above 
the Arctic Circle was completed. 
Owned by the United States Air 
Force, it connects Air Force installa- 
tions in Greenland and Canada. 

The Arctic undersea cable operation 
was carried out under the direction of 


CURRENT 


Topics 
the Western Electric Company with 
the cable placing being handled by the 
Long Lines Department of the 
American Telephone & Telegraph 
Company. Western Electric is 
responsible to the United States Air 
Force for the design, installation and 
testing of rearward communications 
routes for America’s BMEWS project. 

Four nations—Canada, Denmark, 
Great Britain and the United Srates— 
combined to make this first telephone 
cable beneath Arctic waters possible. 

This undersea cable which incor-: 
porates built-in amplifying repeaters 
spaced some 40 miles apart, provides 
reliable, time-tested communication 
channels between the United States 
Air Force Base at Thule, Greenland 
and other military stations in northern 
Canada. It is similar in construction 
to the existing transoceanic telephone 
cables and will be used for military 
purposes. Work on the cable started 
in the summer of 1958 and was com- 
pleted late in 1959 at a cost of about 
$9,000,000. 

It was ‘‘touch and go” as to whether 
the cable would be completed on 
schedule because of two factors: the 
relatively short periods of the year 
when the Baffin Bay and Davis Strait 
region would be free of pack ice and 
icebergs, and the fact that no one had 
ever placed an undersea telephone 
cable in Arctic waters. 

In addition, an emergency devel- 
oped—one with possible serious con- 
sequences for the cable project—in 
June of last year the British cable ship 
Ocean Layer was gutted by fire while 
completing the second transatlantic 
sable. The Ocean Layer had been 
scheduled to lay the second of the 
twin 700-mile BMEWS cables. The 
first cable had been put down success- 
fully the summer before by the British 
Post Office cable ship, HMTS 
Monarch. To surmount the emer- 
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gency, two vessels, the American 
Cable and Radio Corporation's cable 
ship John W. Mackay, and the United 
States Military Sea Transport Service 
cargo ship Lt. James Robinson, worked 
as a team. 

The Mackay, when fully loaded, 
carried only enough cable to reach 
approximately half the distance from 
Cape Dyer to Thule. The other ship 
was used to re-supply her. The 
Robinson was converted from a stand- 
ard cargo ship to a cable carrier by 
around-the-clock shifts of workers at 
an east coast navy shipyard. 

The Mackay loaded cable at the 
Simplex Wire and Cable Company 
docks in Newington, New Hampshire 
in late August then sailed to Cape 
Dyer, on Baffin Island, and began 
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laying cable northward toward Thule. 
Meanwhile, the Robinson loaded the 
remaining lengths of the needed cable 
and sailed directly to Thule. When 
the Mackay completed laying the first 
segment of the cable she buoyed the 
cable end and sailed for Thule. 

There, she reloaded her cable tanks 
from the Robinson, then raced back 
down Baffin Bay until she located, 
picked up, and spliced into the buoyed 
end of the cable. The Mackay then 
continued northwards to Thule, laying 
the last segment of the second cable. 
Several days later she reached there 
and spliced her cable into the shore 
end. 

The first undersea telephone cable 
in Arctic seas had become a reality— 
ahead of schedule. 


YOU CAN ADVANCE SCIENCE EDUCATION 


Today, more than ever before in its 135-year history, there is vital need for 
The Franklin Institute to effectively promote education in science and technology. 
It is imperative that we meet this challenge by providing adequate educational op- 
portunities in these fields. This requires vision, objective planning, and money. 
We have more than enough of the first two requisites, but far too little of the third. 

Our programs are aimed at professional scientists and industry, as well as the lay 
public and young people seeking inspiration and guidance in choosing a career. 
The Institute’s educational programs are impressive, for they begin with students in 
the early grades of our elementary schools and continue throughout an individual's 
professional or industrial life. With more funds at our disposal, the scope and 
vigor of these activities could be greatly increased. 

The Franklin Institute is not richly endowed. 
depending for encouragement and support on an understanding public. 
and conservative management assures wise administration of all funds. 

Your gift or bequest, large or small, will be deeply appreciated and will be used 
effectively to broaden the Institute’s educational usefulness. There is a warm 
satisfaction in giving financial support to an organization that has pioneered in, and 
is dedicated to, the advancement of science and technology through education. 

When property is transferred, title should be in the name of The Franklin 
Institute of the State of Pennsylvania for the Promotion of the Mechanic Arts. 

The Secretary of The Franklin Institute will gladly furnish you with additional 
information. Write to him at The Franklin Institute, Benjamin Franklin Parkway 
at Twentieth Street, Philadelphia 3, Pennsylvania. 


It is a non-profit organization, 
Capable 
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NATIONAL BUREAU OF STANDARDS NEWS 


CONTINUOUS-PRODUCTION VACUUM FURNACE 


A continuous-production vacuum furnace! for evacuating and sealing 
ceramic wafer electron tubes? has been designed by the National Bureau of 
Standards for the Navy Bureau of Aeronautics as part of the Navy’s industrial 
preparedness program. Mechanical development and construction of the fur- 
nace were undertaken on contract by Sanders Associates, Inc., Nashua, N. H. 

With this system, the tube assembly is introduced into a vacuum and 
heated to a high temperature. The combination of heat and vacuum degasses 
the tube, and the high temperature brazes together the metallized surfaces of 
the ceramic wafers to make a hermetically sealed unit. 

The machine was conceived and designed by C. P. Marsden of the Bureau’s 
electron devices laboratory in connection with a study of industrial processes 
suitable for fabricating electronic components.* The automatic production 
methods for ceramic wafer tubes required processing the assemblies into com- 
pletely sealed and ready-to-use units by evacuating, degassing, and sealing in 
a single continuous operation. The vacuum furnace performs these operations 
with control of only temperature, degree of vacuum, and time of processing. 
The high temperature eliminates the need for getters, thus simplifying tube 
design. The furnace has been used to produce up to 200 tube assemblies 
an hour. 

Essentially, the vacuum furnace is a vertical tube, about 11 ft. long, of 
heat-resisting alloy, surrounded by a heating unit, and open to the atmosphere 
at both ends. Through this tube passes a series of hollow cylindrical specimen 
carriers, each separated by a short piston fitting closely inside the tube. These 
parts are introduced into the top of the tube at atmospheric pressure; once 
inside, they are evacuated through ports along the side of the tube and through 
holes in the carriers. The pistons and carriers pass through the tube under 
gravity feed. In the prototype model, the sliding fit of the pistons maintains 
the vacuum at 10-° mm Hg, and the central portion of the tube can be heated 
to 1000° C. After processing, the parts are removed from the lower end of 
the tube at atmospheric pressure and at reduced temperature. 

With the carrier loading and unloading equipment the machine requires 
about 15 ft. vertically. A 15-ft. square of floor area is adequate for the vacuum 
pumps, controls, and other associated equipment. 

The central 7-ft. portion of the tube is surrounded by the electric heating 
element. This element is simply resistance wire supported in the inner side 
of a grooved ceramic cylinder. Five heat reflecting shields surround the 
heating unit; the inner two are of molybdenum, the outer three of stainless steel. 


1 Patent No. 2,894,795, “Automatic Evacuation System for Vacuum Tubes,”’ awarded to 
C. P. Marsden, 14 July ’59. 

2“Ceramic Wafer Tubes for Modular Units,” by C. P. Marsden, Electronics, Vol. 32, 
No. 39 (Sept. 25, 1959). 

§“Project Tinkertoy,”” NBS Tech. News Bul., Vol. 37, p. 11 (1953), 
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The carriers, 2 in. in diameter and about 4 in. long, are of alloy steel and 
fit closely inside the tube. They contain a number of }-in. holes to permit 
escape of gas. The outside bearing surfaces of the pistons are flame-coated 
with tungsten carbide to prevent the bare metal surfaces from brazing together 
in the intense heat of the furnace section. This coating is extremely hard; 
after the machine had been used continuously for 6 months, the coated surfaces 
showed negligible signs of wear. 

Several vacuum pumps are employed to reduce the air pressure near the 
ends of the tube and to maintain a high degree of vacuum in the central heating 
section. Relatively little air leaks into the system past the carriers and 
pistons; the air removed from the central section is mostly that desorbed from 
the surfaces of the parts entering the system. 

A similar machine could well be devised for continuously processing devices 
other than ceramic radio tubes. The coating of lenses or resistors by vapor- 
plating, and the production of vacuum condensers and switches are possible 
applications. Where the processed article must come out of the vacuum near 
room temperature, the lower end of the furnace could be lengthened and 
water-cooled. 


DESIGN OF CORNER REFLECTOR ANTENNAS 


A design analysis of corner-reflector antennas has been carried out by the 
National Bureau of Standards Boulder (Colo.) Laboratories under the sponsor- 
ship of the U. S. Air Force. The results of this evaluation ‘ by H. V. Cottony 
and A. C. Wilson of the Bureau’s Central Radio Propagation Laboratory 
provide comprehensive and detailed information on designing high-perform- 
ance, economical corner reflectors. 

A corner-reflector antenna consists simply of two plane reflecting surfaces 
joined edge to edge to form a corner, which is usually parallel to the ground. 
The driven element is placed in the aperture between the two planes. This 
type of antenna has the advantages of high gain, broad frequency response, 
narrow beam width, low back radiation, low cost, and ease of construction. 

Early experimental work in ionospheric scatter propagation stimulated 
interest in corner reflectors because they are particularly suitable for this mode 
of communication. Even though these antennas had been used prior to the 
new experimental scatter circuits, their performance was known only qualita- 
tively. Therefore as part of its propagation research program, the Bureau 
undertook an investigation of corner deflectors to determine how gain varies 
with changes in width and length of reflecting surfaces, angle of aperture, and 
position of driven element. 

The experimental corner-reflector antenna used for this investigation con- 
sisted of two lattice-type wooden frames, each 12.3 ft. (five wavelengths at 
400 mc.) wide by 12.3 ft. long, supporting the reflecting surfaces. The re- 
flectors were of overlapping strips of sheet aluminum fastened to the frames 
so that removal of one strip would subtract 0.2 wavelength from the surface. 


4“Radiation Pattern of Finite-Size Corner-Reflector Antennas,” by A. C. Wilson and 
H. V. Cottony, TRE Trans. Ant. & Prop. AP-8, No. 2, March 1960; for further technical in- 
formation, see ‘Gains of Finite-Size Corner-Reflector Antennas,” by H. V. Cottony and A. C. 
Wilson, IRE Trans. Ant. & Prop., Vol. AP-6, No. 4, October, 1958. 
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Width could be varied by trimming the lengths of the strips. The two frames 
were pivoted along the same axis, and the angle of aperture could be varied 
from 20 to 180°. 

The driven element was a folded half-wave dipole. Its position could be 
varied from 0.07 to 2.5 wavelengths from the apex of the reflecting surfaces. 
The dipole’s support also served as a balun to transform the impedance of the 
system to about 50 ohm. Fine adjustment, for an exact match to a 50-ohm 
line, was provided by a two-stub tuning unit. Gain was measured at 400 mc 
by comparing the system with a known antenna.°® 

Although in practice collinear arrays of four and more dipoles are used in 
scatter antennas to obtain gains of over 20 db, the general performance can be 
predicted from an experimental system using only a single dipole. Therefore, 
only the single dipole was used for the measurements. 

Results of the measurements indicate that minimum reflector width is 0.5 
wavelength for the smallest usable value of gain. With an increase in width 
up to 2 wavelengths, gain increases. However, beyond this, there is little or 
no increase in gain, and under some circumstances it may actually decrease. 

For maximum gain, the dipole must be in one of several discreet positions— 
referred to as first, second, third, etc., in order of distance from the apex. The 
exact location of each position is a function of the aperture angle. 

When the dipole is in the first position, gain increases monotonically as 
reflecting surface length increases. When the dipole is in the second or third 
positions, the same relation holds true generally, except for surface lengths less 
than 1.5 wavelengths. For these shorter lengths, the relation is not only too 
irregular for simple analysis but also does not provide the antenna designer 
with any useful data. 

Maximum gain along the forward axis seems to be achieved with the dipole 
in the second position. In this case, the reflecting surface should be more than 
two wavelengths long and approximately two wavelengths wide. 

lonospheric scatter antennas designed for the lowest portion of the VHF 
band should have especially low secondary lobes and minimal radiation to the 
rear, as radiation in these undesired directions is likely to result in self-inter- 
ference from multi-path propagation. An important advantage of the corner 
reflector over other types of antennas is that it can be designed to have ex- 
tremely low back radiation. The small amount that does occur results from 
diffraction around the edges and from penetration through small openings in 
the reflector surfaces. 

Radiation patterns have been compiled for a variety of sizes of corner 
reflectors. This material has been further summarized graphically so that 
beam widths and back radiation levels for any combination of surface lengths 
and widths are easily determined. 


5 “Techniques for Accurate Measurement of Antenna Gain,” by H. V. Cottony, NBS 
Circ. No. 598. Available from the Superintendent of Documents, U. S. Government Printing 
Office, Washington 25, D. C. Price, 15 cents. 
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LABORATORIES FOR RESEARCH AND DEVELOPMENT 


announce opportunities for: 


ELECTRICAL ENGINEERS 
SYSTEMS ENGINEERS 
ORGANIC CHEMISTS 

OPERATIONS RESEARCH SPECIALISTS 
MECHANICAL ENGINEERS 
EXPERIMENTAL PSYCHOLOGISTS 
METALLURGISTS 
INDUSTRIAL ENGINEERS 
PHYSICISTS 
AFPLIED MATHEMATICIANS 


to work on challenging problems in research 


Send complete resume to: 
Mr. John E. Christ, Director of Personnel 


THE FRANKLIN INSTITUTE 
Philadelphia 3, Pa. 
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THE FRANKLIN INSTITUTE LABORATORIES 
FOR RESEARCH AND DEVELOPMENT 


In the fields of engineering and the physi- 
cal sciences, a competent and versatile 
staff of several hundred scientists and en- 
gineers, working with modern equipment 
in a creative climate, is trained to bring a 
fresh scientific approach to the solution of 


industrial problems. 
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